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I.         INTRODUCTION 

Combat  modeling  and  simulation  are  important  analysis  tools  available  to  mil- 
itary leaders  and  their  supporting  organizations  which  provide  valuable  insight  into 
the  study  of  military  operations  across  the  spectrum  of  conflict.  Combat  models 
are  employed  frequently  to  conduct  investigative  studies  ranging  from  peacekeeping 
operations,  which  are  relatively  narrow  in  scope,  to  conventional  maneuver  warfare 
executed  at  the  nation-state  level.  Commanders  and  their  supporting  planning  staffs 
regularly  seek  the  resulting  output  from  the  execution  of  combat  models  and  simu- 
lations. Prior  to  the  employment  of  man  and  machine  into  a  hostile  environment, 
tactical  commanders  have  the  capability  to  develop  and  assess  their  maneuver  plans 
by  use  of  combat  models  and  simulations.  Additionally,  combat  models  and  simu- 
lations are  employed  as  training  tools  to  rehearse  command  and  control  actions  and 
operational  procedures.  While  models  and  simulations  are  not  in  themselves  predic- 
tors of  the  outcomes  to  military  operations  or  plans,  their  ability  to  provide  samples 
about  possible  outcomes  is  welcomed  as  a  preliminary  step  prior  to  placing  soldiers 
and  equipment  in  harm's  way.  As  such,  the  designers  of  these  models  constantly 
attempt  to  refine  the  representation  of  warfare  to  ensure  the  model  is  as  accurate 
and  believable  as  possible.  It  is  important  to  note  that  many  of  the  interactions  and 
human  factors  which  are  inherent  to  a  military  force  are  difficult  to  replicate  within 
a  combat  model  or  simulation.  This  fact  must  be  recognized  by  those  employing 
combat  models  and  simulations  for  mission  analysis  and  maneuver  warfare  planning. 

One  of  the  most  important  capabilities  of  combat  models  and  simulations  is 
their  ability  to  calculate  the  attrition  of  personnel  and  systems  of  the  engaged  military 
forces.  The  prediction,  analysis,  and  management  of  the  attrition  to  personnel  and 
equipment  in  combat  are  highly  regarded  as  important  tasks  by  military  leaders  and 
their  staffs.  And,  as  a  result,  the  ability  of  combat  models  and  simulations  to  capture 
attrition  data  is  often  a  measurement  of  the  model's  usefulness.    What  follows  is  a 


discussion  of  the  development  of  attrition  measurements  within  the  context  of  combat 
models  and  simulations. 

A.   BACKGROUND 

Several  important  principles  and  concepts  underlying  the  general  class  descrip- 
tion of  combat  models  and  simulations  are  important  to  this  development  of  attrition 
measurement.  Combat  models  and  simulations  are  generally  subdivided  into  two 
categories,  namely,  high  resolution  and  aggregated. 

High  resolution  models  are  generally  associated  with  smaller  force  compo- 
sitions where  each  combatant  and  its  associated  weapon  or  maneuver  system  are 
explicitly  represented  in  the  simulation.  The  complicated  combat  processes  are  de- 
composed into  highly  detailed  sequences  of  events  and  activities  which  are  carried 
out  explicitly  within  the  simulation.  This  detailed  sequence  of  events  and  activities 
lends  itself  to  simulation  time  steps  which  are  very  short.  It  is  not  uncommon  for 
high  resolution  models  to  represent  time  explicitly  to  the  second,  updating  impor- 
tant information  almost  continuously.  The  goal  of  the  high  resolution  model  is  to 
accurately  model  each  combat  process  so  that  the  simulation  results  are  traceable 
using  mathematical  formulas  and  constructs.  The  characteristics  associated  with  the 
entities  represented  in  the  simulation  are  related  to  their  engineering  specifications, 
physical  performance  data,  and  specified  behavioral  assumptions.  The  existence  of 
this  information  link  from  the  engineering  level  through  the  behavioral  engagement 
assumptions  is  the  single  most  important  quality  of  a  high  resolution  model.  The 
capability  of  representing  weapon  systems  and  equipment  at  the  engineering  level 
enables  an  analyst  to  observe  the  possible  effects  of  alternative  system  specifications 
on  a  pre-defined  simulation.  Individual  system  contributions  to  the  overall  simulation 
are  effectively  analyzed  in  high  resolution  models.  Another  feature  associated  with 
high  resolution  models  is  based  upon  their  reliance  on  stochastic  processes.  Nearly 
all  of  the  simulation  algorithms  above  the  engineering  specification  level  are  proba- 


bilistic  in  nature.  For  example,  consider  a  firing  sequence.  The  firing  system  goes 
through  several  sequential  actions  prior  to  firing  its  weapon.  These  include  target 
acquisition,  target  selection,  firing,  impact,  and  damage  assessment.  As  these  events 
are  not  of  a  fixed  length  in  actual  combat,  their  representation  in  simulation  relies 
upon  the  use  of  highly  probabilistic  random  variables.  Calculations  of  this  nature 
can  lead  to  variance  among  independent  trials  of  a  given  combat  simulation.  In  order 
to  obtain  a  reliable  result  with  respect  to  the  variance  between  stochastic  calcula- 
tions, the  modeler  must  execute  multiple  runs  of  the  same  simulation  to  obtain  the 
"averaged"  result.  One  of  the  most  significant  limitations  of  high  resolution  models 
is  their  sheer  size  and  complexity.  The  computer  codes  supporting  high  resolution 
models  are  highly  complex  and  often  quite  large,  stretching  hardware  and  software  to 
their  limit.  These  computer  codes  are  expensive  to  develop,  run,  and  maintain.  The 
complexity  associated  with  representing  each  individual  entity  of  the  combat  scenario 
at  the  engineering  specification  level  causes  another  significant  limitation  of  high  res- 
olution models,  namely,  force  size.  High  resolution  models  are  limited  to  relatively 
small  forces.  In  ground  maneuver  terms,  this  limitation  is  normally  associated  with 
an  upper  bound  of  roughly  1000  combatants  [Ref.  1]  or  on  the  order  of  battalion  to 
brigade  levels  of  combat. 

Aggregated  models,  as  the  word  implies,  maintain  varying  levels  of  aggregation 
of  the  represented  combatants.  For  example,  an  aggregated  model  may  be  capable 
of  representing  a  group  of  riflemen  with  an  associated  consolidated  firepower  capa- 
bility rather  than  the  individual  soldiers.  The  model  then  portrays  the  aggregated 
capabilities  of  the  group  of  riflemen  represented  graphically  as  a  single  entity.  There- 
fore, an  entity  representing  a  platoon  of  30  riflemen  would  maintain  the  firepower 
and  maneuver  capabilities  associated  with  the  30  individual  weapons  systems  carried 
by  the  30  riflemen.  The  extent  of  aggregation  within  combat  models  varies  based 
upon  the  level  of  combat  being  modeled  and  the  intended  use  of  the  model.  The 
aggregation  process  is  natural  within  the  confines  of  military  combat  models  because 


of  the  hierarchical  command  structure  inherent  to  maneuver  forces.  For  example, 
the  natural  aggregation  process  associated  with  a  generic  ground  maneuver  organiza- 
tion is  individual  combatant,  squad,  platoon,  company,  battalion,  brigade,  division, 
corps,  army,  and  theater.  Depending  upon  the  use  of  the  model  in  question,  one 
of  these  natural  groupings  (normally  no  lower  than  brigade)  would  mark  the  lowest 
level  of  representation  within  the  simulation.  Individual  combatants  and  groups  of 
combatants  which  are  smaller  than  the  basic  size  are  not  explicitly  represented  and 
the  attributes  associated  with  these  combatants  are  lost  to  the  simulation.  Two  im- 
portant and  distinctive  aggregation  patterns  associated  with  models  of  this  type  are 
homogeneous  and  heterogeneous  aggregation,  and,  the  distinction  between  them  is 
particularly  important  for  attrition  calculations.  A  homogeneous  aggregation  is  one 
where  the  combat  power  of  a  unit  with  multiple  and  distinct  combatants  is  combined 
into  a  single  measure  (or  perhaps  one  for  ground  combat  systems  and  a  second  for 
aircraft  systems).  Attrition  computations  are  then  based  upon  the  relative  power  of 
the  aggregated  forces  in  battle  subject  to  their  combat  power  indices.  Contrastingly, 
heterogeneous  aggregation  maintains  a  count  of  the  number  of  surviving  weapons  sys- 
tems of  each  distinct  type.  This  method  of  aggregation  allows  the  model  to  represent 
the  different  levels  of  effectiveness  of  a  particular  weapon  system  against  its  possible 
targets.  This  form  of  aggregation  permits  more  accurate  attrition  modeling  than  the 
homogeneous  model  resulting  in  less  information  loss.  As  a  result,  the  trend  of  most 
modern  aggregated  models  is  to  follow  a  heterogeneous  aggregation  process. 

As  the  size  of  the  modeled  entity  increases  in  aggregated  models,  the  preci- 
sion to  which  they  can  be  located  on  the  simulated  battlefield  may  be  decreased.  A 
modeled  combatant  may  maintain  a  center  location  and  a  frontal  width  where  the 
individual  entities  making  up  the  aggregated  unit  are  assumed  to  be  somewhere  in  the 
unit  region.  Information  about  the  individual  entities  within  the  aggregated  unit  may 
be  kept  in  the  form  of  averages,  but  any  specificity  about  individual  locations  is  lost. 
Thus,  the  first  significant  consequence  of  aggregation  is  the  loss  of  individual  char- 


acteristics  within  aggregated  units  resulting  in  the  use  of  their  averaged  properties 
where  any  inherent  fluctuations  have  been  smoothed.  Analogous  to  the  smoothing 
of  individual  variations  is  the  handling  of  stochastic  processes.  Since  aggregate  mod- 
els do  not  maintain  a  reference  to  what  each  individual  combatant  is  doing  at  any 
given  time,  the  complex  probabilistic  events  found  within  high  resolution  models  are 
replaced  by  averaged  behaviors.  The  firing  process  described  earlier  (target  acquisi- 
tion, target  selection,  firing,  impact,  and  damage  assessment)  within  the  context  of 
high  resolution  models  may  be  replaced  in  an  aggregated  model  by  computing  the 
rate  at  which  an  average  tank  kills  enemy  tanks.  This  derived  attrition  rate  implic- 
itly defines  each  of  the  probabilistic  processes  described  above  with  the  stochastic 
variability  suppressed.  This  fact  leads  to  a  descriptor  generally  associated  with  ag- 
gregated models,  namely,  deterministic.  With  the  engineering  specification  level  of 
detail  and  the  probabilistic  processes  being  suppressed  within  the  aggregated  model, 
many  assumptions  and  scenario  details  get  absorbed  by  the  attrition  rate  number, 
which  in  itself  is  a  highly  variable  entity. 

Lastly,  aggregated  models  lose  the  information  associated  with  event  sequenc- 
ing since  they  do  not  keep  track  of  individual  actions.  Thus,  the  precise  times  of 
critical  events,  such  as  target  kills,  are  not  available  from  aggregated  model  output. 
This  results  in  the  process  averages  being  processed  over  relatively  long  periods  of 
time  (ranging  from  a  few  minutes  to  a  24  hour  period)  where  the  computation  of  total 
casualties  occurs  at  the  end  of  the  time  period.  Again,  careful  consideration  must  be 
given  to  the  attrition  rate  formulation  in  that  there  are  at  least  two  sides  involved  in 
the  exchange  of  fires.  Not  all  combatants  on  either  side  will  survive  to  continue  firing 
for  the  duration  of  the  time  period.  For  example,  consider  the  following  example  and 
resulting  conclusions: 

Extreme  care  must  be  taken  to  make  the  underlying  assumptions  clear 
in  computing  the  rate  data  for  such  models.  For  example,  a  daily  attrition  rate 
for  aircraft  stationed  at  an  airbase  which  is  under  attack  may  be  inconsistent 
with  the  day's  plan  for  those  aircraft  (they  may  be  flying  missions  and  thus 


not  on  the  ground  when  their  base  is  attacked).  A  simulation  which  assesses 
attrition  on  an  aggregated  daily  cycle  cannot  know  whether  the  airbase  attack 
coincided  with  aircraft  on  the  ground,  so  an  assumption  must  be  made  and 
figured  into  the  attrition  rate  data. 

Overall,  then,  the  outstanding  feature  of  aggregation  in  combat  models 
is  information  loss.  As  a  result,  it  is  necessary  to  model  combat  processes  using 
average  behavior  rather  than  individual  behavior.  Although  values  for  the 
average  process  rates  are  generally  computed  outside  the  combat  model,  the 
complex  nature  of  these  rates  demands  that  we  devote  considerable  effort... to 
explaining  their  derivation.  [Ref.  1] 

An  obvious  question  with  regards  to  the  given  strengths  and  limitations  of  the 
described  combat  modeling  and  simulation  approaches  is  why  high  resolution  models 
are  not  extended  to  larger  force  sizes  thus  maintaining  the  admirable  qualities  of  the 
high  resolution  definition.  Unfortunately,  the  difficulty  in  making  this  extension  lies 
within  the  limitations  of  computer  hardware  and  software  and  the  expansive  execu- 
tion time  which  would  be  associated  with  such  large  scale  high  resolution  models. 
The  sheer  number  and  complexity  of  the  individual  combatants  and  their  associated 
characteristics  render  the  use  of  high  resolution  combat  models  for  all  combat  simu- 
lations impractical  and,  in  some  cases,  impossible.  Additionally,  the  complex  nature 
of  command,  control,  and  communications  (C3)  between  echeloned  organizations  is 
extremely  difficult  to  model.  The  synergism  and  inherent  lack  of  measurability  of 
these  traits  would  require  many  assumptions  on  the  part  of  the  modeler  possibly  re- 
ducing the  degree  of  model  believability.  Thus,  the  resulting  segregation  of  combat 
models  into  the  high  resolution  and  aggregated  categories  is  a  matter  of  necessity. 
It  is  important  to  note,  however,  that  the  two  categories  may  not  be  entirely  inde- 
pendent during  their  execution  due  to  the  aggregated  models'  possible  reliance  upon 
the  outcomes  of  high  resolution  battle  simulations  for  the  derivation  of  attrition  rate 
parameters.  While  there  are  several  currently  employed  methods  for  calculating  force 
attrition  within  aggregated  models  (to  be  presented  in  Chapter  III),  it  remains  true 
that  the  attrition  of  personnel,  equipment,  and  supplies  resulting  from  an  engagement 


must,  at  some  point,  be  calculated  within  an  aggregated  simulation.  The  compromise 
which  is  of  specific  interest  to  this  research  is  the  algorithmic  application  of  mathe- 
matical constructs  to  determine  aggregated  force  attrition  based  upon  high  resolution 
results.  For  example,  when  an  aggregated  combat  simulation  reaches  a  point  where 
two  opposing  forces  engage  in  battle,  the  model  may  implicitly  rely  upon  high  res- 
olution simulation  results  generated  under  similar  environmental,  force  composition, 
and  tactical  conditions.  The  end  result  of  the  battle  executed  within  the  aggregated 
simulation  would  be  new  force  levels  for  each  of  the  engaged  combatants  derived  pri- 
marily from  a  deterministic  replay  of  a  representative  battle  using  the  attrition  rate 
coefficients  derived  from  the  high  resolution  simulation  output.  The  intermediate 
process  of  estimating  the  attrition  rates  and  replaying  the  battle  mathematically  is 
the  conceptual  link  between  high  resolution  and  aggregated  models.  With  the  line  di- 
viding high  resolution  and  aggregated  models  clear  in  definition  only,  both  categories 
have  an  associated  collection  of  applications  which,  when  combined,  adequately  meet 
the  needs  of  combat  modeling  and  simulation. 

B.      APPLICABILITY  TO  CURRENT  PROCEDURES 

In  the  United  States  Army  and,  moreover,  the  joint  warfare  community,  one 
of  the  most  frequently  used  high  resolution  models  which  provides  the  data  necessary 
to  generate  the  attrition  rates  for  use  in  aggregated  models  is  called  the  Combat 
Sample  Generator  (COSAGE)  model.  The  COSAGE  model  is  a  (fairly)  high  reso- 
lution, stochastic  model  of  division  level  ground  combat  which  maintains  a  combat- 
ant resolution  to  friendly  platoons  and  enemy  companies.  Within  each  maneuver 
unit,  a  heterogeneous  list  of  individual  weapon  systems  is  arranged  in  combat  forma- 
tions, interactions  between  opposing  weapon  systems  are  computed,  and  individual 
weapons  system  survival  is  governed  by  the  laws  of  probability.  With  respect  to 
aircraft,  close  air  support  of  the  ground  battle  by  tactical  aircraft  is  resolved  down 
to  a  detailed  simulation  of  the  flight  of  individual  aircraft.    At  a  point  in  the  exe- 


cution  of  an  aggregated  model  where  two  opposing  ground  forces  are  positioned  to 
engage  in  combat,  the  representative  battle  simulation  is  carried  out  as  a  mathe- 
matical algorithm  based  primarily  upon  the  results  of  a  high  resolution  simulation 
(COSAGE).  The  killer/victim  results  and  the  associated  time  series  of  individual 
attrition  are  analyzed  and  mathematically  manipulated  resulting  in  the  development 
of  attrition  rate  estimates  for  each  individual  weapon  systems  of  each  engaged  force. 
This  process  is  referred  to  as  the  attrition  calibration  (ATCAL)  Phase  I  process  and 
will  be  discussed  in  detail  in  Chapter  IV.  These  derived  attrition  rates  are  then 
applied  to  a  deterministic  algorithm  which,  when  executed  or  "replayed"  with  the 
force  levels  present  in  the  aggregated  model,  calculates  resulting  force  levels  (by 
system).  These  updated  quantities  are  then  passed  back  to  the  aggregated  model 
as  it  continues  through  its  remaining  processes.  This  secondary  "replay"  process 
is  commonly  referred  to  as  ATCAL  Phase  II.  Currently,  several  theater  level  mod- 
els have  the  capacity  to  employ  ATCAL  Phases  I  and  II.  These  include  Tactical 
Warfare  (TACWAR),  Concepts  Evaluation  Model  (CEM),  and  Integrated  Theater 
Engagement  Model  (ITEM).  Whether  a  particular  theater  level  aggregated  model 
employs  the  ATCAL  algorithm  as  its  method  for  calculating  attrition  is  largely  de- 
pendent upon  the  respective  developers'  and  users'  understanding  of  ATCAL's  theo- 
retical development  and  mathematical  principles,  which,  in  themselves,  remain  to  be 
areas  of  active  research.  See  [Ref.  2]. 

C.      PROBLEM  STATEMENT 

In  order  to  support  current  and  future  investigations  into  the  theoretical  un- 
derpinnings and  mathematical  stability  of  the  defining  algorithms  within  ATCAL 
Phases  I  and  II,  researchers  must  have  a  "user  friendly"  method  of  obtaining  and 
analyzing  the  necessary  high  resolution  output  data  and  validating  the  conjectures 
present  in  said  algorithms.  Additionally,  as  a  starting  point  for  the  mathematical 
analysis,  researchers  must  have  a  thorough  understanding  of  the  basic  principles  as- 


sociated  with  the  various  methods  of  calculating  force  attrition.  This  capability  for 
analysis  and  foundation  for  understanding  do  not  exist  outside  of  the  United  States 
Army  Concepts  Analysis  Agency  (CAA)  and  its  proprietary  implementation  of  the 
COSAGE  model  and  ATCAL  Phase  I  and  II  algorithms.  The  intent  of  this  thesis  is 
to  develop  the  framework  for  a  test  bed  capable  of  conducting  detailed  analysis  of  the 
attrition  problem.  By  developing  realistic  high  resolution  ground  combat  simulations 
using  the  Janus  simulation  software,  capturing  the  detailed  killer/victim  results  and 
the  associated  time  series  of  casualties,  and  applying  these  results  to  the  study  of  the 
aggregated  attrition  problem,  continued  and  specialized  future  research  can  be  con- 
ducted in  this  discipline.  Specifically,  as  an  immediate  result  of  this  research,  future 
inquires  can  be  made  into  the  validity  and  appropriateness  of  the  currently  accepted 
versions  of  the  ATCAL  Phases  I  and  II  algorithms  aiding  in  the  proper  documenta- 
tion of  their  conjectures  and  results.  Moreover,  with  this  vehicle  for  comparative 
analysis,  efforts  can  expand  to  the  development  of  alternative  attrition  algorithms  to 
meet  the  needs  of  future  aggregated  modeling  efforts. 

D.      JANUS  SIMULATION  SOFTWARE 

As  the  initial  step  in  creating  a  functional  test  bed  for  the  study  of  attri- 
tion, one  must  choose  a  vehicle  to  conduct  the  sample  data  base  of  high  resolution 
simulations.  Janus  is  the  primary  high  resolution,  interactive  wargaming  simulation 
software  used  by  the  United  States  Army  for  the  modeling  of  brigade  sized  (and  be- 
low) operations,  and,  its  widespread  availability  within  the  military  community,  ease 
of  use,  and  verified  credibility  make  it  an  obvious  choice  for  meeting  the  demands  of 
this  proposed  analysis.  Named  for  the  two-faced  Roman  god  who  was  the  guardian 
of  portals  and  the  patron  of  beginnings  and  endings,  Janus  had  its  beginnings  as  a 
nuclear  effects  modeling  simulation  developed  by  the  Lawrence  Livermore  National 
Laboratory  of  the  University  of  California.  [Ref.  3]  Its  applicability  to  tactical  train- 
ing and  operations  analysis  quickly  made  it  functional  to  the  United  States  Army 


Training  and  Doctrine  Command  (TRADOC),  Training  Analysis  Center  (TRAC). 
The  current  version  is  intended  to  satisfy  both  the  combat  development  and  training 
communities.  Its  versatility  allows  it  to  be  used  interactively  as  a  man-in-the-loop 
simulation  to  replicate  the  effects  of  command  and  control  on  a  realistic  battlefield 
during  the  simulation  as  enemy  activity  develops.  As  an  analysis  tool,  Janus  can 
be  used  non-interactively  thus  reducing  the  variability  introduced  by  human  decision 
making  during  the  simulation  execution.  Janus  is  a  stochastic  simulation  in  that  the 
system  determines  the  results  of  actions  (like  direct  fire  engagements)  according  to 
laws  of  probability.  While  the  principle  focus  of  Janus  is  on  ground  maneuver  and 
artillery  units,  it  also  models  weather  and  its  effects,  day  and  night  visibility,  engineer 
support,  minefield  employment  and  breaching,  rotary  and  fixed  wing  aircraft,  resup- 
ply  operations,  and  chemical  environments.  The  inherent  flexibility  of  this  simulation 
software  allows  its  user  to  accurately  model  a  multitude  of  combat  scenarios. 

E.       SCOPE  AND  LIMITATIONS 

This  research  is  unique  in  that  there  has  been  no  other  attempt  to  establish 
an  independent  capability  for  analyzing  the  aggregated  attrition  problem,  specifically 
the  ATCAL  Phase  I  and  II  algorithms.  This  work  will  validate  the  capability  of  con- 
structing a  high  resolution  ground  combat  simulation  in  Janus  for  the  purpose  of 
deriving  attrition  rate  parameters  using  the  ATCAL  Phase  I  and  II  algorithms.  Ad- 
ditionally, there  has  been  no  attempt  to  configure  the  Janus  output  data  to  conform 
to  the  requirements  of  the  current  ATCAL  Phase  I  and  II  algorithms  (currently  exist- 
ing in  the  FORTRAN  computer  language).  This  procedure  and  the  transformation  of 
the  output  data  is  practical  only  within  the  confines  of  analysis  and  algorithm  valida- 
tion because  of  the  time  expenditure  required  to  develop  the  ground  combat  scenario 
in  Janus  and  the  nonexistence  of  automated  software  to  link  the  Janus  killer/victim 
results  to  current  aggregated  models.  With  initial  efforts  confined  to  the  manual  com- 
pilation of  Janus  output  data  and  the  explicit  formulation  and  calculation  of  basic 
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mathematical  models,  this  research  could  be  expanded  to  focus  on  the  compatibility 
and  automation  of  the  data  linkage  between  Janus  the  current  and/or  future  ATCAL 
computer  codes.  However,  this  is  not  the  focus  of  this  thesis.  And,  as  only  one  specific 
ground  combat  simulation  will  be  presented  as  part  of  this  thesis,  follow-on  work  is 
applicable  in  the  development  of  additional  scenarios  with  varying  mission  objectives, 
force  definitions,  environmental  conditions,  and  terrain  selections  in  order  to  general- 
ize results  otherwise  limited  to  specific  simulation  parameters.  A  robust  collection  of 
simulations  would  be  required  to  properly  validate  the  conjectures  presented  in  the 
problem  statement.  As  such,  the  results  presented  in  this  work  relate  only  to  the 
specific  scenario  and  situations  modeled  and  should  not  be  considered  to  be  the  final 
solution. 


11 


12 


II.         THE  SIMULATION  SCENARIO 

In  order  to  generate  the  required  high  resolution  output  data  with  Janus  for 
the  development  of  attrition  rate  numbers  using  the  ATCAL  Phase  I  and  II  algo- 
rithms (analogous  to  the  process  described  in  Chapter  I,  Section  B),  the  author  has 
chosen  to  simulate  a  heavy  brigade  combat  team  (BCT)  combat  scenario.  The  typical 
heavy  BCT  includes  a  heterogeneous  mix  of  direct  fire,  indirect  fire,  and  supporting 
maneuver  systems.  Generically,  the  force  would  include  tanks,  armored  personnel  car- 
riers, artillery  systems,  air  defense  artillery  systems,  engineering  equipment,  wheeled 
vehicles,  air  assets,  and  obstacles.  The  BCT  is  the  preferred  level  of  combatant  or- 
ganization for  this  study  because  of  the  described  combined  arms  composition  and 
the  ability  to  portray  this  force  in  Janus.  As  a  high  resolution  model,  Janus  can 
effectively  support  a  force  type  of  this  size  and  capture  the  detailed  results  of  all 
direct  and  indirect  fire  engagements.  With  these  results,  the  attrition  rate  estimation 
process  (ATCAL  Phase  I)  and  the  deterministic  battle  replay  (ATCAL  Phase  II)  can 
be  conducted. 

The  tactical  scenario  chosen  for  this  combat  simulation  is  one  of  practical 
relevance.  As  this  research  focuses  primarily  on  the  modeling  of  ground  combat 
as  part  of  an  overlying  larger  tactical  operation,  the  author  has  chosen  to  simulate 
the  mechanized  task  force  operations  in  desert  environmental  conditions  as  a  base- 
line scenario.  As  the  Janus  simulation  package  can  be  an  interactive  training  tool 
when  employed  as  a  man-in-the-loop  simulation,  it  is  often  used  by  TRADOC  at  the 
United  States  Army's  combat  arms  training  centers.  Specifically,  the  Armor  Officer 
Advanced  Course  (AOAC)  located  at  the  Unites  States  Army  Armor  Center,  Fort 
Knox,  Kentucky,  employs  Janus  as  a  tool  to  train  and  evaluate  the  mission  planning 
and  tactical  decision  making  skills  of  its  officer  students.  With  a  specified  mission  of 
training  the  current  and  future  leaders  of  the  United  States  maneuver  warfare  orga- 
nizations, AOAC  and  its  focus  on  modern  ground  maneuver  tactics,  techniques,  and 
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procedures  provide  an  excellent  simulation  opportunity  for  this  research.  The  Janus 
simulation  used  during  this  research  was  executed  by  AOAC  Class  97-3  at  the  Fort 
Knox  Brigade  Battle  Simulation  Center.  [Ref.  4]  And,  the  simulation  is  available 
for  review  through  coordination  with  Professor  Bard  K.  Mansager,  Department  of 
Mathematics,  Naval  Postgraduate  School.  The  stochastic  nature  of  the  simulation  is 
defined  not  only  by  the  probabilistic  events  which  occur  as  a  function  of  the  intrinsic 
algorithms,  but  also  by  the  human  factors  associated  with  the  particular  officer  stu- 
dents who  executed  the  simulation.  If  multiple  runs  of  this  scenario  were  executed, 
each  would  produce  different  results  based  upon  these  stochastic  processes. 

A.      MECHANICS  OF  THE  SIMULATION 

At  the  highest  level  of  organization,  Janus  allows  the  definition  of  six  combat- 
ant sides  where  each  side  could  represents  an  autonomous  force.  Under  each  side, 
forces  can  be  further  organized  into  task  forces  where  modeled  weapons  and  maneu- 
ver systems  can  be  organized  and  manipulated  in  groupings  that  mirror  their  real  life 
task  organization.  After  the  defined  combatants  are  organized  into  task  forces,  one  or 
more  of  these  task  forces  are  output  to  a  specific  computer  workstation  for  user  con- 
trol. Janus  permits  up  to  24  computer  workstations  to  be  in  use  during  a  simulation. 
The  24  workstations  are  divided  into  side  1  through  side  6  terminals  based  upon  the 
number  of  users  assigned  to  manipulate  the  respective  sides  and  task  forces.  [Ref.  5] 
The  simulation  used  for  this  research  makes  use  of  sides  1  and  2  and  will  hereafter 
be  referred  to  as  BLUE  and  RED,  respectively.  This  assignment  represents  the  roles 
of  the  combatants  during  the  simulation.  In  this  simulation,  BLUE  represents  the 
friendly  Unites  States  forces  being  controlled  by  the  AOAC  97-3  student  officers  and 
RED  represents  a  Soviet-style  enemy  forces  being  controlled  by  designated  opposing 
force  personnel. 
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1.  BLUE  Force  Composition  and  Organization 

The  BLUE  force  (BLUEFOR)  is  made  up  of  the  United  States  maneuver 
and  weapons  systems  appearing  in  Appendix  A,  Table  XIII  on  page  64;  and,  when 
organized  for  combat,  these  systems  represents  a  heavy  BCT.  The  BLUEFOR  is 
considered  to  be  "balanced"  in  terms  of  its  armor  to  mechanized  infantry  ratio.  As 
an  integrated  combined  arms  fighting  force,  the  BLUEFOR  systems  are  organized  into 
a  typical  BCT  hierarchical  command  and  control  structure  as  detailed  by  Table  I.  For 
the  purposes  of  the  AOAC  student  officers  who  executed  the  Janus  simulation,  this 
fictitious  force  was  designated  as  the  1st  Brigade  Combat  Team  (1  BCT). 

Unit  Description  Major  Systems 

4  Direct  Fire  Maneuver  Battalions  M1A1,  M2 

2  Indirect  Fire  Artillery  Battalions  M109A6 

1  Air  Defense  Artillery  Company  BSFV 

1  Chemical  Support  Company  IF  TRU 

1  Engineer  Support  Battalion  AVLM 

1  Logistics  Support  Battalion  IF  TRU 

1  Tactical  Air  Control  Party  A10  (6) 

Table  I.  BLUEFOR  Task  Organization  for  AOAC  97-3  Janus  Simulation 
[Ref.  6] 

2.  RED  Force  Composition  and  Organization 

The  RED  force  (REDFOR)  is  made  up  of  Soviet  style  equipment  listed  in 
Appendix  A,  Table  XII  on  page  63;  and,  when  organized  for  combat,  it  represents 
a  less  than  full  strength  Soviet  style  motorized  rifle  regiment  [MRR(-)].  Note  the 
greater  than  three  to  one  ratio  of  mechanized  infantry  to  armor,  implying  that  the 
REDFOR  is  better  equipped  for  mechanized  infantry  operations  than  tank  warfare. 
Additionally,  the  REDFOR  has  a  significant  number  of  indirect  fire  artillery  systems 
at  its  disposal.  The  defensive  REDFOR  MRR(-)  is  organized  within  the  hierarchical 
command  and  control  structure  detailed  by  Table  II.  For  the  purposes  of  the  AOAC 
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student  officers  who  executed  the  Janus  simulation,  this  fictitious  enemy  force  was 
designated  as  the  168th  MRR(-). 

Unit  Description  Major  Systems 

2  Direct  Fire  Maneuver  Battalions  BMP 

1  Direct  Fire  Tank  Battalion  T80 

4  Indirect  Fire  Artillery  Battalions  2S1,  2S5,  BM-21 

1  Air  Defense  Artillery  Company  SA7,  ZSU-23 

1  Engineer  Support  Company  MTK-LI 

1  Logistics  Support  Company  IF  TRU 

1  Helicopter  Company  HIND  (6) 

Table  II.  REDFOR  Task  Organization  for  AOAC  97-3  Janus  Simulation 

[Ref.  6] 


B.      TERRAIN  REPRESENTATION 

The  Janus  simulation  uses  digitized  terrain  developed  by  the  United  States 
National  Imaging  and  Mapping  Agency  (NIMA).  The  terrain  selected  for  this  sim- 
ulation is  located  within  the  United  States  Army's  National  Training  Center  (NTC) 
near  Barstow,  California,  covering  a  100  km2  area  with  actual  military  grid  coordi- 
nate locations  defined  by  the  corner  points  NK800400,  NK790400,  NK790390,  and 
NK800390.  The  terrain  is  displayed  in  Janus  similar  to  that  used  by  military  person- 
nel during  actual  training  and  combat  operations.  It  includes  contour  lines  (elevation 
information),  roads,  rivers,  vegetation,  and  urban  areas.  The  visibility  issues  as- 
sociated with  terrain  elevation  and  visual  obstructions  are  realistically  portrayed  in 
Janus.  These  modeled  attributes  will  degrade  (or  enhance)  a  combatant's  line-of-sight 
(LOS)  and  movement.  The  modeled  vehicles  and  equipment  must  maneuver  under 
the  conditions  represented  by  the  digital  terrain  and  its  overlying  natural  and  man- 
made  features.  For  example,  only  those  vehicles  with  a  "swim"  capability  for  crossing 
water  obstacles  are  permitted  to  do  so  without  employing  bridging  assets;  and,  when 
executed,  the  vehicles  will  proceed  at  a  reduced  movement  rate. 
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The  environmental  conditions  and  time  of  battle  are  also  controlled  in  Janus. 
The  user  can  specify  environmental  conditions  including  temperature,  cloud  condi- 
tions, and  air  characteristics  affecting  smoke  propagation.  Weather  and  its  effects 
influence  the  visibility  and  movement  of  the  modeled  combatants.  The  start  time 
can  be  set  in  order  to  execute  the  combat  simulation  at  a  particular  time  of  day. 
This  accounts  for  simulating  day  and  night  operations  which  also  affect  visibility  and 
movement  attributes.  The  environmental  conditions  used  during  this  simulation  were 
summer,  desert  with  14  kilometers  (km)  visibility,  and  a  0800  start  time.  Since  Janus 
accurately  models  the  BLUEFOR  and  REDFOR  weapons  systems  as  a  function  of 
each  system's  capabilities  as  defined  in  the  database,  users  must  consider  all  military 
aspects  of  employing  their  forces  just  as  they  would  in  actual  combat.  Invariably, 
Janus  will  reveal  the  weaknesses  associated  with  poor  planning  and  reward  those  who 
have  tactically  sound  plans.  As  the  combat  scenario  for  this  research  was  executed 
by  students  of  the  mounted  warfare  profession,  the  author  assumes  the  employment 
of  tactically  sound  techniques,  procedures,  and  plans. 

C.      SCENARIO  EXECUTION  SUMMARY 

The  initial  dispositions  of  BLUEFOR  and  REDFOR  are  depicted  in  Ap- 
pendix A,  Figures  1  and  2  on  pages  65  and  67,  respectively.  The  BLUEFOR  mission 
is  to  attack  in  zone  to  destroy  or  cause  the  withdrawal  of  the  168th  MRR(-)  secu- 
rity zone  forces  who  are  controlling  the  key  terrain  in  and  around  grid  coordinate 
NK620130  (Avawatz  Mountain  Passes)  [Ref.  6].  The  BLUEFOR  executes  a  sup- 
porting attack  in  the  north  followed  by  a  main  attack  in  the  south.  The  REDFOR 
conducts  a  defensive  operation  with  an  objective  of  stopping  the  penetration  of  the 
BLUEFOR  and  defending  the  significant  terrain  as  mentioned  above  [Ref.  6].  The 
REDFOR  executes  a  layered  defense  with  two  distinct  defensive  belts  supported  by 
deliberate  minefields,  obstacles,  and  barriers. 
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At  the  onset,  the  BLUEFOR  is  positioned  to  conduct  offensive  operations 
using  two  distinct  avenues  of  approach  (A A).  For  the  purpose  of  this  research,  the 
northern  and  southern  AAs  will  be  named  AA  North  and  AA  South,  respectively.  The 
main  effort  of  the  BLUEFOR  attacks  along  A  A  South  with  the  supporting  attack 
along  AA  North.  As  a  matter  of  mounted  warfare  tactics,  the  BLUEFOR  leads 
with  reconnaissance  assets  in  order  to  locate  and  direct  fires  upon  the  REDFOR's 
defensive  assets.  Following  closely  behind  are  the  BLUEFOR's  initial  maneuver  and 
engineering  assets.  This  battlefield  organization  will  support  direct  fire  engagements 
with  the  REDFOR  upon  their  detection  by  the  BLUEFOR  reconnaissance  elements. 
The  BLUEFOR  engineering  elements  are  positioned  forward  to  facilitate  the  clearing 
of  obstacles  and  breaching  of  REDFOR  minefields.  Next  in  the  BLUEFOR  order  of 
battle  is  the  bulk  of  the  fighting  forces  primarily  consisting  of  tanks  and  armored 
personnel  carriers.  The  indirect  fire  artillery  assets  follow  in  the  BLUEFOR  combat 
formation.  They  are  positioned  in  an  attempt  to  be  within  range  of  the  REDFOR 
while  maintaining  reasonable  security  from  REDFOR  counter-artillery  operations. 
Finally,  the  support  oriented  vehicles  and  equipment  make  up  the  rear  of  the  attacking 
BLUEFOR.  Positioned  well  to  the  rear  and  in  support  of  the  BLUEFOR,  six  A- 10 
close  air  support  (CAS)  aircraft  are  available  to  be  employed  as  a  coordinated  part 
of  the  attack. 

The  REDFOR,  while  maintaining  a  layered  defense  supported  by  maneuver 
obstacles  and  minefields,  also  maintains  a  counterattack  capability  with  reserve  forces 
located  to  the  rear  of  their  main  defensive  area.  These  combatants  can  be  employed 
to  reinforce  weakened  defensive  areas  or  to  counter  a  penetrating  attack  made  by  the 
BLUEFOR.  As  detailed  by  Table  II,  the  REDFOR  is  primarily  a  mechanized  force 
consisting  of  armored  personnel  carriers.  As  such,  the  REDFOR's  tanks  are  employed 
judiciously  and  in  positions  where  they  can  have  the  most  significant  effect  on  the  at- 
tacking BLUEFOR.  One  of  the  most  significant  capabilities  of  the  REDFOR  is  found 
in  the  number  of  supporting  indirect  fire  artillery  assets.    The  REDFOR  maintains 
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both  a  close  and  long  range  indirect  fire  capability  while  keeping  the  majority  of  its 
artillery  assets  outside  of  the  maximum  range  of  the  BLUEFOR's  artillery  units.  The 
REDFOR  emphasizes  the  use  of  its  indirect  fire  systems  and  their  associated  radar 
systems  to  locate  and  conduct  counter-battery  fire  against  the  BLUEFOR.  As  the 
simulation  progresses,  the  REDFOR  attempts  to  halt  the  advance  of  the  attacking 
BLUEFOR. 

1.       Execution  Summary:  0  to  30  Minutes 

The  BLUEFOR  begins  offensive  operations  by  employing  its  leading  recon- 
naissance assets  in  AA  North.  Shortly  thereafter,  the  BLUEFOR  and  REDFOR 
make  contact  with  each  other  and  exchange  direct  fires.  At  approximately  11  min- 
utes, the  REDFOR  long  range  artillery  assets  fire  special  munitions  placing  two  area 
minefields  each  containing  480  anti-armor  mines  to  thwart  the  advance  of  the  attack- 
ing BLUEFOR  in  the  north.  As  of  yet,  the  BLUEFOR  attacking  along  AA  South 
have  not  made  contact  with  the  REDFOR  defensive  positions  in  that  area.  During 
and  immediately  following  the  BLUEFOR's  minefield  breaching  operations  in  the 
north,  a  heated  direct  fire  battle  takes  place.  Both  sides  employ  their  indirect  fire 
artillery  assets  as  their  reconnaissance  assets  fix  each  other's  locations  on  the  battle- 
field. With  28  BLUEFOR  and  9  REDFOR  combatants  killed  during  this  15  minute 
localized  battle,  the  BLUEFOR  attacking  along  A  A  North  are  able  to  overwhelm  the 
REDFOR  defensive  forces  and  break  through  the  first  defensive  line  in  the  north. 
Shortly  thereafter,  at  approximately  27  minutes,  the  REDFOR  long  range  artillery 
assets  reposition  to  alternate  firing  positions  for  security  purposes. 

See  Appendix  A,  Figure  3  on  page  69  for  a  pictorial  representation  of  the 
battlefield  at  the  beginning  of  the  scenario.  Table  III  below  details  the  force  attrition 
at  the  completion  of  the  0  to  30  minute  time  period.  See  Appendix  A,  Figure  4  on 
page  71  for  a  pictorial  representation  of  the  battlefield  at  30  minutes. 
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BLUEFOR  Losses  REDFOR  Losses 
Direct  Fire                                                                   34  10 

Indirect  Fire  7  4 

Minefields  4  0 

Cumulative  Losses  at  30  Minutes  45  14 

Table  III.  Force  Losses  for  AOAC  97-3  Janus  Simulation  at  30  Minutes 

2.        Execution  Summary:  31  to  60  Minutes 

The  BLUEFOR,  who  met  with  success  in  the  north,  proceed  east-northeast 
on  a  high  speed  A  A  and  are  met  with  little  resistance  from  the  REDFOR.  At  ap- 
proximately 40  minutes,  the  attacking  BLUEFOR  moving  along  AA  South  come  in 
contact  with  the  first  layer  of  REDFOR  defenders  in  the  south.  The  majority  of 
the  BLUEFOR  following  the  two  attacks  begin  movement  north-northeast  to  take 
advantage  of  the  successful  penetration  of  the  REDFOR's  northern  defenses.  The 
BLUEFOR  in  contact  in  the  south  encounter  and  begin  breaching  operations  on 
three  linear  minefields  positioned  in  front  of  the  REDFOR.  As  BLUEFOR  recon- 
naissance forces  identify  REDFOR  targets  to  their  front  they  direct  indirect  artillery 
fires  against  them  resulting  in  only  marginal  success.  At  approximately  45  minutes, 
the  remainder  of  the  REDFOR's  defending  forces  in  the  first  defensive  layer  engage 
the  attacking  BLUEFOR  from  AA  South.  The  REDFOR  defense  in  the  south  is 
much  more  formidable  than  that  which  was  found  in  the  north.  The  BLUEFOR  has 
some  difficulty  identifying  penetration  points  through  the  REDFOR's  defenses.  The 
BLUEFOR  suffers  a  large  number  of  casualties  during  this  battle.  The  main  body  of 
the  BLUEFOR  compresses  into  a  more  compact  movement  formation  and  increases 
its  movement  rate. 

Table  IV  details  the  force  attrition  at  the  completion  of  the  31  to  60  minute 
time  period.  See  Appendix  A,  Figure  5  on  page  73  for  a  pictorial  representation  of 
the  battlefield  at  60  minutes. 
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BLUEFOR  Losses  REDFOR  Losses 
Direct  Fire                                                                   51  14 

Indirect  Fire  9  4 

Minefields  4  0 

Cumulative  Losses  at  60  Minutes  64  18 

Table  IV.  Force  Losses  for  AOAC  97-3  Janus  Simulation  at  60  Minutes 

3.        Execution  Summary:  61  to  90  Minutes 

By  65  minutes,  the  BLUEFOR  attacking  in  AA  South  are  massed  at  the  for- 
ward edge  of  the  linear  minefields  and  conducting  passage  through  the  cleared  lanes. 
At  the  same  time,  six  A-10  CAS  aircraft  enter  an  aerial  holding  area  for  possible 
employment  against  the  REDFOR  in  the  first  defensive  belt.  The  REDFOR  has  a 
small  ambush  or  counter-attack  force  made  up  of  two  tanks  and  two  armored  person- 
nel carriers  positioned  at  the  extreme  south-southeast  of  its  first  defensive  layer  in 
order  to  block  the  advancing  BLUEFOR  from  using  a  high  speed  AA.  However,  the 
overwhelming  BLUEFOR  combat  power  moving  through  this  area  quickly  engages 
and  destroys  this  team.  At  this  point,  the  BLUEFOR  begins  to  reorient  its  forward 
movement  to  the  north  to  take  advantage  of  the  previous  success  there.  At  approxi- 
mately 83  minutes,  the  BLUEFOR  attacking  to  the  north-northeast  along  AA  North 
make  contact  with  the  second  layer  of  REDFOR  defenses.  The  BLUEFOR  begins  an 
indirect  artillery  barrage  against  these  REDFOR. 

Table  V  details  the  force  attrition  at  the  completion  of  the  61  to  90  minute 

time  period.    See  Appendix  A,  Figure  6  on  page  75  for  a  pictorial  representation  of 

the  battlefield  at  90  minutes. 

BLUEFOR  Losses  REDFOR  Losses 
Direct  Fire                                                                   65  25 

Indirect  Fire  10  5 

Minefields  7  0 

Cumulative  Losses  at  90  Minutes  82  30 

Table  V.  Force  Losses  for  AOAC  97-3  Janus  Simulation  at  90  Minutes 
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4.  Execution  Summary:  91  to  120  Minutes 

At  approximately  93  minutes,  the  six  A-10  CAS  aircraft  depart  the  aerial  hold- 
ing area  to  attack  the  second  layer  of  the  RED  FOR  defenses  from  the  southeast  to 
the  northwest.  Shortly  thereafter,  at  100  minutes,  two  REDFOR  radar  controlled 
anti-aircraft  systems  engage  and  destroy  three  A-10  CAS  aircraft.  The  remaining 
aircraft  adjust  their  course  and  continue  their  attack  on  the  defending  REDFOR 
combatants.  Prior  to  their  egress,  the  remaining  three  A-10  CAS  aircraft  are  de- 
stroyed by  REDFOR  shoulder-fired  anti-aircraft  missiles  located  in  the  northeast.  By 
105  minutes,  the  main  BLUEFOR  attacking  elements  clear  the  linear  minefields  and 
being  their  movement  north-northeast  along  the  high  speed  AA.  At  1 16  minutes,  these 
attacking  BLUEFOR  make  contact  and  exchange  initial  fires  with  the  remainder  of 
the  REDFOR's  second  layer  defensive  forces  in  the  central  and  southern  portion  of 
the  maneuver  area.  Heavy  exchanges  of  artillery  occur  during  the  development  of 
this  high  intensity  battle. 

Table  VI  details  the  force  attrition  at  the  completion  of  the  91  to  120  minute 
time  period.  See  Appendix  A,  Figure  7  on  page  77  for  a  pictorial  representation  of 
the  battlefield  at  120  minutes. 

BLUEFOR  Losses  REDFOR  Losses 
Direct  Fire                                                                  92  35 

Indirect  Fire  24  9 

Minefields  12  0 

Cumulative  Losses  at  120  Minutes  128  44 

Table  VI.  Force  Losses  for  AOAC  97-3  Janus  Simulation  at  120  Minutes 

5.  Execution  Summary:   121  to  150  Minutes 

During  the  battle  within  the  second  layer  of  the  REDFOR  defense,  six  ad- 
ditional A-10  CAS  aircraft  are  called  forward  from  their  staging  area  to  attack  the 
defending  forces.  The  majority  of  the  six  aircraft  are  highly  successful  until  their 
destruction  at  131  minutes  by  radar  controlled  anti-aircraft  fires.  Heavy  exchanges  of 
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direct  fires  continue  to  occur  and  the  attacking  BLUEFOR  divides  its  formation  into 
two  distinct  attacks,  one  to  the  north-northwest  and  another  along  the  high  speed 
A  A  to  the  east-northeast. 

Table  VII  details  the  force  attrition  at  the  completion  of  the  121  to  150  minute 
time  period.  See  Appendix  A,  Figure  8  on  page  79  for  a  pictorial  representation  of 
the  battlefield  at  150  minutes. 

BLUEFOR  Losses  REDFOR  Losses 
Direct  Fire                                                                    139  49 

Indirect  Fire  32  9 

Minefields  22  0 

Cumulative  Losses  at  150  Minutes  193  58 

Table  VII.  Force  Losses  for  AOAC  97-3  Janus  Simulation  at  150  Minutes 

6.        Execution  Summary:   151  to  180  Minutes 

The  BLUEFOR  attack  to  the  north-northwest  meets  with  a  tank  heavy  REDFOR 
defensive  element  including  some  reserve  forces  coming  from  behind  the  second  de- 
fensive layer.  The  BLUEFOR  reacts  by  passing  additional  tank  assets  through  to 
the  front  of  its  formation  in  an  attempt  to  counter  the  REDFOR  tank  threat  with 
like  firepower.  At  158  minutes,  the  BLUEFOR  in  the  northwest  directs  indirect  ar- 
tillery fires  against  the  REDFOR  with  little  effect.  At  165  minutes,  the  main  body 
of  the  BLUEFOR  which  was  following  the  lead  attacking  forces  moves  farther  for- 
ward to  concentrate  fires.  The  direct  fire  battle  in  the  north-northwest  becomes 
a  battle  of  attrition  with  the  BLUEFOR  maintaining  numerical  superiority.  How- 
ever, this  translated  into  a  target  rich  environment  for  the  REDFOR.  A  successful 
penetration  of  the  defense  occurs  and  BLUEFOR  continues  their  northwestern  move- 
ment past  the  REDFOR's  second  layer  of  defense.  As  the  breakthrough  occurs,  the 
REDFOR  in  the  south  take  advantage  of  the  flanking  shot  opportunities  against  the 
passing  BLUEFOR  and  are  rewarded  with  a  large  number  of  kills.  This  time  period 
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is  marked  by  the  BLUEFOR's  largest  number  of  direct  fire  losses  as  compared  to  any 
other  30  minute  time  period. 

Table  VIII  details  the  force  attrition  at  the  completion  of  the  151  to  180  minute 
time  period.  See  Appendix  A,  Figure  9  on  page  81  for  a  pictorial  representation  of 
the  battlefield  at  180  minutes. 

BLUEFOR  Losses  REDFOR  Losses 
Direct  Fire                                                                  254  63 

Indirect  Fire  44  9 

Minefields  27  0 

Cumulative  Losses  at  180  Minutes  325  72 

Table  VIII.  Force  Losses  for  AOAC  97-3  Janus  Simulation  at  180  Minutes 

7.       Execution  Summary:  181  to  211  Minutes 

The  BLUEFOR  divides  into  two  relatively  equally  sized  formations  with  one 
attacking  north-northeast  and  the  second,  south  of  the  first,  attacking  along  the  high 
speed  A  A  to  the  northeast.  At  188  minutes,  the  REDFOR  commits  its  reserves  in 
the  north  to  counter  the  BLUEFOR  attack  in  that  sector.  After  some  success,  the 
numerically  superior  BLUEFOR  destroys  the  employed  reserves.  At  198  minutes,  the 
BLUEFOR  in  the  south  completely  clears  the  second  defensive  layer  and  all  of  its 
associated  minefields.  At  206  minutes,  the  REDFOR  commits  addition  reserve  forces 
against  the  BLUEFOR  attacking  in  the  south  along  the  high  speed  AA  in  an  attempt 
to  thwart  the  advance.  At  211  minutes,  the  combat  simulation  is  terminated  simulat- 
ing the  achievement  of  one  or  more  of  the  BLUEFOR  commander's  ending  conditions. 
Battle  termination  standards  may  include  attrition  to  a  specified  force  level,  attrition 
to  a  specified  number  of  remaining  maneuver  systems,  and/or  the  accomplishment  of 
stated  objectives.  This  particular  combat  simulation  was  undoubtedly  stopped  due  to 
a  self-imposed  time  constraint  as  the  simulation  was  executed  by  and  for  the  benefit 
of  AOAC  Class  97-3. 
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Table  IX  details  the  force  attrition  at  the  completion  of  the  simulation.  See 
Appendix  A,  Figure  10  on  page  83  for  a  pictorial  representation  of  the  battlefield  at 
the  completion  of  the  simulation. 

BLUEFOR  Losses  REDFOR  Losses 
Direct  Fire                                                                                     270  79 

Indirect  Fire  45  9 

Minefields  29  0 

Cumulative  Losses  at  Completion  Minutes  344  88 

Table  IX.  Force  Losses  for  AOAC  97-3  Janus  Simulation  at  Completion 


D.      UTILIZATION  OF  THE  SIMULATION  OUTPUT 

One  of  the  most  advantageous  properties  of  the  Janus  simulation  software  is 
its  ability  to  capture  information  about  the  execution  of  the  combat  scenario.  The 
force  loss  totals  detailed  in  the  tables  above  are  just  samples  of  the  information 
available  to  the  Janus  analyst.  The  Janus  post  execution  output  files  which  are 
important  to  the  attrition  rate  analysis  of  this  research  include  the  direct  fire  report, 
coroner's  report,  killer/victim  scoreboard,  and  the  artillery  impacts  report.  While 
several  other  reports  are  available,  these  four  reports  will  provide  the  time  series  of 
casualties  data  required  for  the  development  of  attrition  rate  estimates  applicable  to 
the  ATCAL  Phase  I  algorithm.  With  the  required  simulation  output  data  in  hand, 
what  remains  to  be  developed  are  the  mathematical  constructs  which  serve  as  the 
basis  for  the  ATCAL  Phase  I  and  II  algorithms. 
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III.         ATTRITION  METHODOLOGIES  FOR 
AGGREGATED  MODELS 

A.   BACKGROUND 

Prior  to  any  detailed  analysis  of  the  attrition  rate  estimate  process  as  it  applies 
to  the  ATCAL  Phase  I  and  II  algorithms  (the  primary  goal  of  this  research),  it  is 
important  to  establish  a  foundation  with  respect  to  the  various  methods  for  calculat- 
ing attrition  and  determining  the  outcomes  of  battles  within  aggregated  simulations. 
In  general,  nearly  all  attrition  computation  methodologies  can  be  categorized  by  two 
techniques:  the  Anderson  (Force-Ratio  [FR])  model  and  some  type  of  Lanchester- 
type  model.  Integral  to  the  second  attrition  approach  are  the  Bonder-Ferrell  method 
and  the  Combat  Analysis  (COMAN)  maximum  likelihood  estimation  (MLE)  method 
for  determining  attrition  coefficients  for  use  in  Lanchester-type  equations  (of  which 
variants  are  used  in  the  ATCAL  algorithms).  Although  the  COMAN  MLE  ap- 
proach is  not  specifically  used  within  any  current  models,  its  applicability  to  the 
attrition  modeling  cannot  be  overstated.  Each  of  these  distinct  attrition  methodolo- 
gies maintains  a  relevance  to  modern  aggregated  modeling.  For  example  (although 
not  all-inclusive),  the  Anderson  (FR)  method  is  one  of  the  two  options  for  assessing 
attrition  in  TACWAR  4.1,  the  Bonder-Ferrell  method  is  employed  within  the  United 
States  Army's  Vector  In  Commander  (VIC)  model,  and  the  ATCAL  methodology  is 
used  in  CEM  and  ITEM.  (TACWAR  also  has  an  option  to  use  the  ATCAL  attrition 
methodology  rather  than  FR.)  One  of  the  most  interesting  facts  of  this  research  is 
that  ATCAL  uses  none  of  the  attrition  methodologies  (explicitly)  as  they  are  de- 
scribed below.  Rather,  ATCAL  uses  an  ad  hoc  method  for  attrition  rate  estimation. 
This  is  truly  one  of  the  mysteries  of  ATCAL  for  which  this  research  attempts  to 
establish  the  vehicle  for  further  analysis. 
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B.      THE  ANDERSON  (FR)  BASED  ATTRITION  MODEL 

In  this  theory,  a  heterogeneous  force  (its  overall  strength)  can  be  described 
by  one  scalar  value  representing  the  consolidated  individual  systems'  contributions 
to  the  force  as  a  whole.  Each  system  maintains  an  assigned  value  (score)  relative 
to  some  common  weighting  criteria.  Summing  over  all  systems  multiplied  by  their 
respective  value  (score)  results  in  the  generation  of  an  overall  firepower  index  (FPI) 
for  the  force  being  represented.  As  this  FPI  is  a  single  scalar  value,  the  resulting 
process  is  characterized  as  a  homogeneous  force  ratio  attrition  model.  The  explicit 
representation  of  the  heterogeneous  force  is  replaced  by  this  single  scalar  measure 
of  combat  power.  This  strategy  can  be  expressed  mathematically  for  a  hypothetical 
force,  X,  as  shown  below  in  Equation  III.  1 . 


FPIX  =  YlSfXi  (Firepower  Score)  (III.l) 

i  =  integer  index  each  representing  a  unique  weapon  system  type,  total  of  m 

5t-  =  value  (score)  of  one  weapon  system  of  type  i 

Xi  =  number  of  weapon  system  of  type  i 

Supposing  that  the  FPI  of  two  opposing  forces  (X  and  Y)  are  represented  by 

FPIX 
FPlx  and  FPIy,  respectively,  then  the  force  ratio,  FR^y  =  ,  gives  a  measure  of 

FPIy 

relative  combat  power  at  the  instant  at  which  the  values  are  calculated.  The  FR^-y 
may  be  used  to  determine  force  missions,  unit  postures,  casualty  computations,  and 
mission  success.  Specifically,  the  FR^y  can  be  indexed  against  personnel  loss  curves 
to  determine  attrition  to  modeled  forces.  The  loss  curves  employed  in  this  process 
are  generally  derived  from  the  results  of  historical  battles.  This  fact  is  considered  to 
be  one  the  inherent  strengths  of  the  Anderson  (FR)  method.  On  the  other  hand,  one 
of  the  perceived  weaknesses  of  this  method  is  associated  with  the  determination  of 
the  individual  weapon  system  values  (scores)  as  discussed  below. 
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1.        Determination  of  Weapon  System  Value 

There  is  no  general  agreement  in  the  modeling  community  on  how  to  deter- 
mine the  individual  weapon  system  values  (scores),  Sf.  The  fundamental  methods 
for  determining  firepower  score  values  include  analysts'  perception  of  combat  value 
(the  naive  estimate),  historical  combat  performance,  measures  of  weapon  firepower 
and  support  platform  characteristics,  and,  finally,  measures  of  what  a  weapon  can 
kill.  This  last  method  attempts  to  compute  the  value  (score)  for  a  weapon  as  being 
proportional  to  the  rate  at  which  it  destroys  the  value  of  opposing  enemy  systems, 
and,  as  such,  becomes  a  circularly  defined  problem.  Moreover,  this  approach's  de- 
pendence on  the  "rate"  at  which  weapons  destroy  opposing  weapons  in  battle  infers 
the  reliance  upon  assumed  or  derived  non-negative  scalar  attrition  rates  (one  for  each 
weapon  system  firing  at  every  opposing  weapon  system  for  both  forces)  which  account 
for  the  highly  detailed  engagement  process  described  in  Chapter  I.  The  mathematics 
of  this  method,  called  the  AntiPotential-Potential  (APP)  method,  amounts  to  solv- 
ing for  the  individual  weapon  system  scores  through  a  linear  system  of  equations  (in 
matrix  form)  as  detailed  below  by  Equations  III. 2  and  III. 3.  The  solutions  to  these 
equations  are  found  by  solving  an  equivalent  system  of  eigenvalue  problems  [Ref.  1], 
for  which  the  development  is  beyond  the  scope  of  this  discussion.  However,  as  a  prac- 
tical example,  Appendix  B  details  the  application  of  the  Anderson  (FR)  method  and 
the  subordinate  APP  process  to  hypothetical  high  resolution  output  data  representa- 
tive of  the  AO  AC  97-3  Janus  simulation  described  in  Chapter  II.  The  calculations  are 
sampled  from  the  heterogeneous  battle  between  selected  BLUEFOR  and  REDFOR 
weapon  systems.  Note  that  this  computational  excursion  does  not  account  for  every 
weapon  system  modeled  in  the  Janus  Simulation  in  an  attempt  to  keep  the  associated 
eigenvalue  problem  numerically  solvable  in  closed  form  without  the  aid  of  automated 
software. 
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CXS?    =    EK^Sj  (III.2) 

m 


i=l 


Force  X  =  {xi,  £2,  #3, . . . ,  xm}  m  distinct  weapon  systems 

Force  Y  =  {yi,  t/2,  V3,  ■  •  ■ ,  yn}  n  distinct  weapon  systems 

S*  =  value  of  one  type  i  system  in  X  force 

Sj  =  value  of  one  type  j  system  in  Y  force 

Ktj  =  matrix  of  rates  at  which  one  Xi  system  kills  ^systems 

Lji  =  matrix  of  rates  at  which  one  Yj  system  kills  ^systems 

Cx  =  proportionality  constant 

Cy  =  proportionality  constant 

2.       Advantages  and  Disadvantages 

The  use  of  historical  data,  both  in  the  generation  of  personnel  and  equipment 
loss  curves  and  in  the  determination  of  weapon  system  values,  generates  a  certain  level 
of  model  believability  and  personal  comfortableness  with  the  Anderson  (FR)  method 
of  attrition.  Moreover,  aggregated  simulations  using  this  method  are  credited  with 
reproducing  comparable  results  to  the  recorded  outcomes  of  actual  battles.  This 
trait  is  often  considered  important  when  attempting  to  validate  and  verify  a  model's 
usefulness.  The  desire  to  extend  these  result  to  the  modeling  of  future  battles  is  a 
task  which  generates  a  considerable  amount  of  discussion.  With  respect  to  determin- 
ing weapon  system  values  (scores),  the  APP  method  and  its  associated  eigenvalue 
problem  are  an  attractive  procedure  for  calculating  homogeneous  force  values  and 
resulting  force  ratios  for  use  in  aggregated  models  because  of  their  mathematical  ease 
of  computation  (regardless  of  their  appropriateness). 
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Conversely,  there  are  a  significant  number  of  limitations  inherent  to  this  pro- 
cedure. A  major  criticism  of  this  method  is  that  the  APP  scoring  method  may  be 
inconsistent  with  the  personnel  loss  curves  which  are  generally  based  upon  historical 
data.  Additionally,  since  Equation  III.  1  is  additive  across  weapon  systems,  there  is 
no  ability  to  represent  the  interactions  among  complementary  weapons  on  the  battle- 
field. This  is  an  especially  serious  shortcoming  considering  the  nature  of  the  current 
battlefield  and  the  synergies  of  modern  weapon  systems.  One  must  exercise  caution 
when  conducting  comparative  assessments  of  force  mix  or  force  balance  since  the 
scalars  FPI*  and  FPIy  do  not  discriminate  between  balanced  and  unbalanced  forces 
(in  terms  of  size).  An  FPI^  =  100  and  an  FPIy  =  100  are  equal  forces  using  this 
method  regardless  of  the  number  of  systems  represented  or  their  respective  values. 
The  reader  will  also  note  that  Equation  III.  1  is  linear  in  the  number  of  weapons  X{ 
of  each  weapon  type  i.  Therefore,  a  model  employing  the  Anderson  (FR)  method 
cannot  be  expected  to  effectively  answer  specific  questions  about  the  contributions  of 
individual  weapons.  For  example,  the  ramifications  of  an  unmanageably  large  num- 
ber of  a  specific  weapon  may  be  intuitively  obvious,  but,  the  FR  model  will  linearly 
increase  the  overall  FPI  without  regard  to  practicality.  It  is  the  opinion  of  the  author 
that  the  overwhelmingly  mysterious  characteristic  of  the  Anderson  (FR)  method  lies 
within  the  development  of  the  individual  weapon  system  values  (scores).  The  reader 
may  have  surmised  that  in  order  to  apply  the  Anderson  (FR)  method,  the  analyst  is 
forced  to  arbitrarily  assume  values  (scores)  for  each  of  the  modeled  weapon  systems 
or  to  rely  upon  the  killer/victim  results  of  high  resolution  simulation  to  develop  the 
attrition  rates  necessary  to  employ  the  APP  method  of  solving  for  weapon  system 
values  (scores).  Associated  with  using  the  high  resolution  simulation  results  to  derive 
the  kill  rates  within  the  APP  method  for  increased  believability,  the  reader  should  be 
aware  of  some  of  the  inherent  weaknesses  which  plague  the  APP  eigenvalue  problem. 
The  score  values  derived  from  the  eigenvalue  problem  are  extremely  dependent  upon 
scenario  characteristics  and  the  number  of  weapons  systems  which  translate  into  tar- 
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get  engagement  opportunities.  Changes  to  the  attrition  rates  (the  key  elements  to 
the  APP  method)  will  typically  cause  all  weapon  system  score  values  to  change  for 
both  forces  in  ways  that  are  hard  to  describe  or  justify.  Consider  another  author's 
evaluation  of  the  APP  method  of  determining  weapon  system  values: 

The  eigenvector  score  values  sometimes  change  in  ways  that  are  hard 
to  explain  and  that  have  been  called  paradoxical  by  some.  For  example,  a 
shift  in  fire  distribution  designed  to  increase  the  kill  rates  of  a  higher  value 
enemy  target  can  sometimes  reduce  the  total  value  index  of  the  firing  force. 
Whether  this  is,  in  fact,  paradoxical  depends  on  how  deeply  the  relationships 
are  followed.  If  there  are  few  of  the  high  value  target,  then  shifting  fire  away 
from  a  lower  value  but  more  numerous  enemy  system  might  very  well  result 
in  a  lower  total  value  being  killed  by  the  firing  force. 

Other  anomalies,  however,  are  harder  to  explain.  The  numeric  values 
of  the  scores  are  sometimes  oversensitive  to  small  changes  in  the  input  kill  rate 
matrices...  [Ref.  1] 

It  is  with  this  knowledge  that  the  analyst  is  drawn  to  developing  similar 
methodologies  using  nonlinear  value  equations,  a  formidable  excursion  which  is  be- 
yond the  scope  of  this  research. 

3.       Typical   Implementation   of  the   Anderson    (FR) 
Method 

A  generic  aggregated  combat  model,  like  TACWAR  may  employ  the  Ander- 
son (FR)  method  for  attrition  calculations  in  the  following  way.  Assume  that  a  hypo- 
thetical model  represents  opposing  forces  down  to  the  division  level  of  organization. 
The  model  would  compute  attrition  for  the  divisions  in  a  combat  sector  with  an  up- 
date cycle  equal  to  the  overall  time  cycle  of  the  model,  possible  every  12  or  24  hours. 
Prior  to  the  calculation  of  attrition,  the  simulation  schedules  several  events  to  update 
the  forces.  These  tasks  may  include  movement  of  forces  into  and  out  of  the  combat 
sectors,  allocation  of  reinforcements,  determination  of  which  force  is  the  attacker  for 
the  cycle,  determination  of  the  defender's  posture  for  the  cycle,  and  computation  of 
the  results  of  any  scheduled  air  battles  for  the  cycle.  At  a  point  during  the  execu- 
tion of  a  generic  aggregated  combat  simulation  where  the  conditions  are  met  for  two 
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opposing  divisions  to  engage  in  battle,  the  Anderson  (FR)  method  would  be  invoked 
to  determine  the  respective  forces'  FPIs.  In  a  simplified  sense,  the  resultant  FR^y 
determines  the  outcome  of  the  battle  with  respect  to  losses.  After  computing  (mostly 
estimating)  the  parameters  required  for  the  Anderson  (FR)  and  subordinate  APP 
calculations,  the  resulting  FRj^y  values  are  indexed  against  the  appropriate  person- 
nel casualty  curve  to  determine  overall  attrition  levels.  (Unfortunately,  indexing  the 
FR^y  against  the  casualty  rate  and  equipment  loss  curves  may  not  be  consistent 
due  to  differing  computational  methods.  The  FR^y  is  derived  using  a  linear  model 
including  the  APP  process,  while  the  casualty  rate  and  equipment  loss  curves  may 
have  been  derived  nonlinearly  from  historical  data.)  During  each  time  cycle  of  the 
simulation  (possibly  12  or  24  hours),  the  APP  eigenvalue  problem  must  be  recursively 
solved  to  determine  the  updated  weapon  system  scores  and  the  resulting  forces'  FPIs 
due  to  their  dependence  on  the  number  of  individual  weapon  systems  remaining.  This 
entire  process  is  simply  repeated  as  the  simulation  continues  through  additional  time 
cycles. 

C.      LANCHESTER-TYPE  DIFFERENTIAL  ATTRITION 
MODELS 

Prior  to  discussing  the  Bonder-Ferrell  and  COMAN  MLE  methods  for  esti- 
mating attrition  rates,  the  reader  must  gain  a  basic  understanding  of  the  overlying 
framework  to  which  the  attrition  rates  are  applied.  Lan Chester- type  models  origi- 
nated with  the  early  works  of  F.W.  Lanchester,  a  British  engineer  and  inventor,  in 
an  attempt  to  quantify  the  advantages  of  concentrating  fires  in  modern  warfare.  He 
hypothesized  that  ancient  warfare  was  merely  a  collection  of  one-on-one  duels  where 
force  size  primarily  determined  battle  outcomes.  The  resulting  mathematical  model 
representing  this  ancient  form  of  warfare  is  referred  to  as  Lanchester's  Linear  Law 
or  Firer-Target  |  Firer-Target  (FT|FT)  attrition.  On  the  other  hand,  Lanchester 
supposed  that  modern  warfare  maintained  the  ability  to  concentrate  the  effects  of 
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multiple  weapons  on  surviving  targets,  thus  creating  a  many-on-one  condition.  In  a 
hypothetical  battle  between  forces  X  and  Y,  this  condition  would  make  the  force  X 
casualty  rate  proportional  to  the  number  of  surviving  firing  systems  in  force  Y.  And, 
the  resulting  casualty-exchange  ratio  would  be  inversely  proportional  to  the  force 
ratio.  This  mathematical  model  for  modern  warfare  conditions  is  termed  Lanch- 
ester's  Square  Law  or  Firer  |  Firer  (F|F)  attrition.  These  underlying  Lanchester 
models,  with  various  extensions  and  improvements,  supply  the  basic  framework  to 
the  ATCAL  Phase  I  and  II  algorithms,  which  will  be  discussed  in  Chapter  IV. 

1.       Lanchester's  Linear  Law  (FT|FT  Attrition) 

The  primary  assumption  of  ancient  warfare  is  that  any  weapon  system  (soldier) 
not  engaged  in  combat  must  wait  until  an  enemy  soldier  becomes  available  before 
joining  the  battle.  The  assumption  implies  that  the  number  of  soldiers  of  force  X  put 
forth  on  the  battlefield  to  engage  force  Y  is  determined  by  the  number  of  soldiers  of 
force  Y  present.  Therefore,  under  the  conditions  of  ancient  warfare,  there  should  be  no 
advantage  to  concentrating  forces.  A  particular  soldier  from  a  force  is  either  engaged 
in  a  one-on-one  duel,  or  not.  Intuitively,  each  force  size  is  a  function  of  time,  where 
the  continuous  real  time  variables  x(t),  y(t),  and  t  are  approximations  to  their  discrete 
analogs  in  an  actual  battle.  The  change  of  each  force  size  with  respect  to  time  depends 
upon  the  number  of  firers  engaged  in  combat  and  can  be  represented  mathematically 

by  the  (homogeneous)  Lanchester  Linear  Law  (FT|FT),  Equations  III. 4  and  III. 5, 

dx  dy 

where  —  and  —  represent  the  change  in  the  size  of  forces  X  and  Y,  respectively, 

over  time. 

dx 

—    =    -ax(t)y(t)  with  x(0)  =  x0  (III.4) 

dt 

J    =    -bx(t)y(t)  with  2/(0)  =  y0  (III.5) 

x(t)    =    size  of  force  X  at  time  t 
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y(t)    =    size  of  force  Y  at  time  t 

a    =    attrition  rate  at  which  an  individual  Y  firer 

kills  X  targets 
6    =     attrition  rate  at  which  an  individual  X  firer 

kills  Y  targets 

2.  Lanchester's  Square  Law  (F|F  Attrition) 
Extending  the  assumptions  of  ancient  warfare  to  modern  battles,  one  gains  the 
attribute  that  many  firers  can  engage  a  single  target,  thus  achieving  the  advantages 
associated  with  concentration  of  fires.  Two  general  cases  of  combat  developed  out 
of  this  study,  namely,  area  fire  and  aimed  fire.  First,  Lan Chester  assumed  area  fire 
to  be  a  situation  where  a  force  spreads  its  fires  over  a  general  area  occupied  by  the 
enemy  rather  that  aiming  directly  at  individual  enemy  weapon  systems.  Assuming 
that  these  area  fires  are  uniformly  distributed  over  the  target  area  and  that  the 
area  is  independent  of  the  number  of  targets,  the  Lanchester  Linear  Law  (FT|FT) 
was  assumed  to  hold  true  for  area  fires.  Second,  aimed  fire  assumes  that  individual 
targets  are  identified  and  attacked  by  one  or  more  opposing  firers.  In  this  case, 
Lanchester  assumed  that  the  attrition  of  force  X  depended  on  how  many  force  Y 
firers  were  directing  aimed  fires  at  force  X.  Mathematically,  the  basic  homogeneous 


form  of  Lanchester's  Square  Law  (F|F)  can  be  described  by  a  system  of  differential 

dx 
equations,  Equations  III. 6  and  III. 7,  where  — 

at 

size  of  forces  X  and  Y,  respectively,  over  time. 


dx  dv 

equations,  Equations  III. 6  and  III. 7,  where  —  and  —  represent  the  change  in  the 

dt  dt 


dx 

—    =     -ay{t)  with  x(0)  =  x0  (III.6) 

J     =     -bx{t)  w\thy(0)  =  yo  (HI.7) 


x{t)    =    size  of  force  X  at  time  t 
y(t)    =    size  of  force  Y  at  time  t 
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a    =    attrition  rate  at  which  an  individual  Y  firer 

kills  X  targets 
b    =    attrition  rate  at  which  an  individual  X  firer 

kills  Y  targets 

The  homogeneous  representation  of  Lanchester's  Square  Law  (F|F)  can  be  ex- 
tended to  meet  the  needs  of  heterogeneous  forces  by  accounting  for  all  of  the  possible 
firer/target  combinations  in  the  simulated  battle.  The  resulting  heterogeneous  for- 
mulation of  Lanchester's  Square  Law  (F|F)  is  provided  by  the  system  of  differential 
equations,  Equations  III. 8  and  III. 9,  shown  below. 


dr-  n 

-i  =  -£M^(t)                                            (IH.8) 

ai  j=\ 

dv-  m 

=  -£***;.•*.■(*)                                            (HI.9) 


13 

dt 


2  =  1 


X{(t)  =  size  of  force  X,  system  i,  at  time  t 

yj(t)  =  size  of  force  Y,  system  j,  at  time  t 

dij  —  attrition  rate  at  which  Y3  firers  kill  X{  targets 

iftij  =  fraction  of  Yj  firers  allocated  against  Xi  targets 

bji  =  attrition  rate  at  which  Xi  firers  kill  Yj  targets 

tyji  =  fraction  of  Xi  firers  allocated  against  Y3  targets 

X  force  =  {xi,X2,X3,...,xm} 

Y  force  =  {yi,tfe,]fey. .  .,yn} 

The  reader  should  be  aware  that  Equations  III. 8  and  III. 9  represent  one  pos- 
sible form  of  the  Lanchester  Square  Law  (F|F)  equations.  While  the  equations  pre- 
sented represent  work  which  was  primarily  accomplished  during  the  1960s,  the  overall 
methodology  has  evolved  over  time.  In  the  1970s,  fire  allocation  factors  which  influ- 
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ence  the  attrition  rates  were  developed  by  submodels  executed  in  support  of  the 
overall  process.  Nevertheless,  the  equations  described  above  appear  in  a  fundamental 
form  which  aids  in  the  understanding  of  this  topic. 

3.        The  Significance  of  aZJ  and  bji 

As  Lanchester's  original  purpose  was  a  qualitative  justification  of  the  principle 
of  concentration  as  an  advantageous  tactic  in  modern  warfare,  the  derivation  of  the 
heterogeneous  attrition  rate  coefficients  a^  and  bji  was  of  no  particular  concern.  How- 
ever, today,  as  Lanchester  based  attrition  models  are  employed  as  part  of  theater  level 
aggregated  models,  the  significance  of  these  attrition  rates  cannot  be  overemphasized. 
As  part  of  the  model  verification  process,  current  aggregated  models  must  provide 
reasonable  and  quantitative  results.  These  numerical  values  assigned  to  atJ  and  67, 
are  what  represent  the  synergies  of  battle  between  various  weapons  systems  and  the 
very  nature  of  specific  actions  on  the  battlefield.  Their  derivation  is  what  brings 
believability  to  the  equations  present  within  combat  models.  Several  methods  have 
been  proposed  and  utilized  over  the  years  to  derive  the  attrition  coefficients.  The 
remaining  background  discussion  of  this  chapter  focuses  on  two  of  these  methods  of 
attrition  rate  estimation  which  are  currently  in  use. 

D.      THE  BONDER-FERRELL  METHOD  FOR  ATTRITION 
RATE  ESTIMATION 

The  Bonder-Ferrell  method  for  estimating  attrition  rates  makes  use  of  an  in- 
dependent analytical  model  which  does  not  depend  on  input  from  other  sources. 
Instead,  a  stochastic  process  model  is  developed  for  each  Yj  firing  at  a  single  passive 
type  Xi  target,  and  the  expected  time  between  casualties  (£[T^])  is  determined,  where 
T-j  is  a  random  variable  representing  the  time  it  takes  a  single  Yj  firing  system  to 
kill  one  Xi  target.  Assuming  that  the  casualty  process  is  a  Markov  (Poisson)  process 

(see  [Ref.  7]),  the  attrition  rate  parameters  can  then  be  expressed  as  atj  =  and 

E[Tij\ 
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bji  =  .  While  Bonder's  original  work  focused  on  deriving  the  equation  to  com- 

pute  E[T]  by  calculating  the  distribution  of  the  number  of  shots  required  to  achieve 
a  specified  number  of  hits,  Ferrell's  extension  of  this  work  modeled  this  process  as  a 
Markov  renewal  process  (see  [Ref.  7])  with  three  discrete  states  (new  engagement,  hit, 
and  miss)  and  the  associated  transition  probabilities  and  times.  Thus,  the  expected 
time  required  to  kill  an  enemy  target  is  the  long  run  expected  time  required  within 
the  Markov  renewal  process  until  moving  on  to  a  new  engagement.  As  early  studies 
in  this  area  focussed  on  homogeneous  forces,  continued  work  in  the  1970s  extended 
the  concept  to  heterogeneous  battles  by  combining  conditional  kill  rates  with  target 
priority  lists.  Additionally,  target  acquisition  was  characterised  as  either  parallel  or 
serial. 

One  of  the  most  significant  advantages  to  the  Bonder-Ferrell  method  of  at- 
trition coefficient  estimation  is  that  the  computations  are  based  upon  what  is  going 
on  in  the  battle  at  the  given  moment.  There  is  no  requirement  to  reference  a  data 
base  library  of  high  resolution  simulation  output  to  match  up  the  conditions  of  battle 
between  the  overlying  model  and  the  high  resolution  battle.  Additionally,  since  the 
underlying  data  used  in  these  computations  are  basically  measurable,  the  results  are 
more  credible  to  the  user.  This  method  is  generally  considered  to  be  restrictive  in 
that  the  in-depth  engagement  model  is  analytic,  thus  suppressing  some  detail.  This 
form  of  attrition  rate  estimation  also  precludes  the  ability  to  capture  synergistic  ef- 
fects of  weapon  system  combinations  as  the  underlying  process  does  not  account  for 
weapon  system  interactions.  The  reliance  upon  engineering  level  data  does  lend  this 
approach  to  the  explicit  definitions  of  engagements;  however,  the  amount  of  data 
required  on  each  weapon  system  is  great.  Relying  upon  closed  form  equations  for  the 
development  of  the  atj  and  6JZ ,  the  Bonder-Ferrell  method  allows  for  the  re-evaluation 
of  the  attrition  coefficients  at  each  time  step  of  the  simulation  by  simply  recalling  and 
recomputing  the  required  input  values. 
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E.  THE  COMAN  MAXIMUM  LIKELIHOOD  ESTIMA- 
TION (MLE)  APPROACH  TO  ATTRITION  RATE 
ESTIMATION 

The  COMAN  MLE  attrition  estimation  methodology  is  a  fitted  parameter 
model  which  examines  the  time  series  of  casualty  times  (from  a  high  resolution  sim- 
ulation) and  computes  the  maximum  likelihood  estimate  for  the  mean  time  between 
casualties.  While  other  statistical  methods  exist  for  determining  point  estimates  of 
this  type  (including  the  method  of  moments,  Bayesian  estimation,  and  the  method  of 
least  squares),  the  MLE  approach  is  the  only  method  which  has  had  any  significant 
application  in  combat  simulations  [Ref.  8].  Similar  to  the  Bonder- Ferrell  approach, 
the  COMAN  MLE  method  relies  upon  the  assumption  that  the  time  between  casu- 
alties is  a  Markov  (Poisson)  process  with  the  memoryless  property  (time  until  the 
next  casualty  does  not  depend  upon  the  previous  time  before  the  last  casualty).  The 
COMAN  MLE  method  is  most  easily  explained  by  detailing  its  mathematical  con- 
struct. The  underlying  assumption  of  this  method  is  that  the  high  resolution  simula- 
tion data  represents  a  sample  from  the  continuous  analogue  to  the  deterministic  homo- 
geneous equations  of  Lanchester's  Square  Law  (F|F)  (see  Equations  III. 6  and  III. 7), 
i.e.,  the  continuous-parameter  Markov  chain  with  transition  probabilities  given  by 
Equations  III.  10  and  III.  11  [Ref.  9]. 

Prob[X  casualty  in  a  small  time  interval  of  length  At]    =    an(t)At 

(111.10) 

Prob[K  casualty  in  a  small  time  interval  of  length  At]    =    bm(t)At 

(III.ll) 

n(i)    =    the  number  of  Y  combatants  at  time  t 
m(t)    =    the  number  of  X  combatants  at  time  t 
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Additionally,  consider  the  following  random  variable  definitions  and  their  as- 
sociated properties  required  for  the  development  of  the  MLE  attrition  coefficient  es- 
timation process.  The  reader  can  assume  that  any  realization  of  a  described  random 
variables  will  be  notationally  represented  by  its  corresponding  lower-case  letter. 


c?   =   i 


cl 


1  if  the  kth  casualty  is  an  X  combatant, 

0  otherwise, 

1  if  the  k     casualty  is  an  Y  combatant, 
0  otherwise, 


K 
Cj     =    ^2  ck  total  number  of  X  casualties 
k=i 

K 

c?    =    2J  ck  total  number  of  Y  casualties 
fc=i 

Cj  +  Cj    =    K  total  number  of  X  and  Y  casualties 

nk      the  realization  of  the  number  of  Y  combatants 
n(tk)     = 

just  after  the  occurrence  of  the  k     casualty 

rrik      the  realization  of  the  number  of  X  combatants 
m(tk)     = 

just  after  the  occurrence  of  the  kth  casualty 

Interpreting  the  random  variable  definitions  given  above,  note  that  during 
the  time  interval  [tk,tk+i]  there  are  m,k  X  combatants  and  nk  Y  combatants  alive 
for  values  of  k  =  (0, 1,2, . . . ,  K  —  1).  Armed  with  high  resolution  simulation  output 
data  of  the  form  {(ii,  Cj  ,  cx ),  (<2,  Cj  ,  c2  ), . . . ,  (£#,  cK,  cK)},  the  statistical  estimates 
for  the  continuous  time  Markov  chain  representation  of  the  deterministic  form  of 
Lanchester's  Square  Law  (F|F)  can  be  computed  using  the  method  of  maximum 
likelihood  estimation  (MLE).  As  a  battle  develops,  the  casualty  processes  for  the  X 
and  Y  forces  are  simply  two  superimposed  Poisson  processes  with  respect  to  time  (a 
random  variable).  Developing  and  maximizing  these  likelihood  functions  from  their 
applicable  joint  probability  density  function  (PDF)  (representing  the  times  between 
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casualties)  will  result  in  the  derivation  of  the  attrition  rate  estimates,  a  and  6.  The 
maximum  likelihood  estimation  process  can  be  summarized  by  the  steps  detailed  in 
Table  X. 


Step  1 


Find  the  probability  density  function  for  the  time  until  an  X  casualty. 
Find  the  probability  density  function  for  the  time  until  a  V  casualty. 


Step  2 


Construct  the  respective  likelihood  functions. 

(i.e.,  the  density  functions  for  the  observed  sequences  of  events) 


Step  3 


Determine  the  values  of  a  and  b  that  maximize  the 
respective  likelihood  functions. 


[Ref.  9] 


Table  X.  Maximum  Likelihood  Estimation  Process  Steps 


1.       Maximum  Likelihood  Estimation  Steps 

Applying  the  three  maximum  likelihood  estimation  steps  to  the  stochastic  rep- 
resentation of  Lanchester's  Square  Law  (F|F)  represented  by  Equations  III. 10  and  III. 11 
will  result  in  the  statistical  derivation  of  a  and  b.  The  author  has  chosen  to  detail 
the  development  of  the  estimated  attrition  parameters  a  and  6  from  the  X  casualty 
point  of  view.  Initially,  C*  and  1  —  C*  will  be  used  to  represent  X  and  Y  casualty 
occurrences,  respectively. 

a.  Step  1 

Let  the  random  variable  Sx  denote  the  time  between  X  casualties  with 
the  well-known  exponential  PDF  fsx(sx)  =  an  exp~^an^Sx .  Similarly,  the  PDF  for 
the  time  between  Y  casualties  is  fsY{sy)  —  bm  exp~^bm^Sy .  Further,  let  the  random 
variable  S  denote  the  time  between  any  two  consecutive  casualties  be  defined  as  the 
minimum  of  Sx  and  Sy-  The  resulting  PDF  describing  the  time  between  any  two 
casualties  is 

fs(s)  =  (an  +  6m)  exp-(an+6m)s .  (111.12) 
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Incorporating  the  previously  defined  variable  C*,  which  equals  one  when  the  kth 
casualty  is  an  X  combatant  and  zero  when  the  kih  casualty  is  a  Y  combatant,  one 
can  further  specify  Equation  III.  12  above.    Note  that  C*  =  1  if  Sx  <  Sy.    The 

/777 

probability  of  this  occurrence  is  Prob(Cfc    =  1)  = — .   Likewise,  a  Y  casualty 

an  +  bm 

hrrt 

occurs  if  Sx  >  Sy  which  implies  C*  =  0  and  Prob(Cjf  =  0)  = — .  Therefore, 

an  +  bm 

the  joint  probability  of  the  random  variables  S  and  C*  is 


Prob(S<«,Cjf  =  z)    =    (1  -  exp-(an+6m)s)  •  Prob(cf  =  i) 

where  0  <  s  <  oo  and  i  =  0, 1 


.  l—SL-yi  ■  (-^— )'-<?  •  (i  -  exp-(°»+6""s). 

an  -\-  bm  an  -\-  bm 


Therefore,  the  representative  joint  PDF  of  S  and  C*  is 


'«*<*•*>  =  (^r^-(^^-(-+6m)^"+M< 


=    (an)^(6m)1-c?exp-(an+6m)s  (111.13) 

where  c*  =  0, 1. 


b.  Step  2 

The  likelihood  function  is  developed  by  considering  that  the  K  casu- 
alties have  occurred  at  times  (<i,^2^3i  •  •  •  ^k)  with  a  total  number  of  X  and  Y  ca- 
sualties represented  by  c?£  and  c^,  respectively,  where  Cj  +  c^  =  K.  Recalling  that 
Cl?  =  1  «=»  C\  -  0  and  Cjf  =  0  <*=>  Cf  =  1,  Equation  III.  13  can  be  functionally 
written  as 

(anfc_1)c?(&m*_1)(*  eXp-(an*-1+6mfc-l)(*k-'*-l) .  (111.14) 
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Equation  III.  14  represents  the  contribution  of  the  kth  casualty  to  the  likelihood  func- 
tion. As  the  times  between  casualties  are  independent  random  variables,  the  likeli- 
hood function  for  the  observed  sequence  of  events,  Equation  III.  15,  is  the  product  of 
all  of  the  independent  contributions. 


K 

n 


L{a,b)  ](a^_1)^(6mjt_1)c^exp-(an*-1+bmfc-l)(<':-tfc-l) 

(111.15) 


c.  Step  3 

Finally,  the  COMAN  MLE  process  requires  the  maximization  of  this 
likelihood  function.  However,  it  is  a  common  mathematical  practice  to  maximize  the 
natural  logarithm  (In)  of  the  likelihood  function  as  detailed  below  by  Equation  III. 16. 
Both  maximum  values  will  occur  at  the  same  point  [Ref.  10]  and  the  logarithmic  form 
is  more  easily  manipulated. 


K  k  K 

ln[L(a,6)]    =    Yl  ck  ln(an*_i)  +  J2  cl  ln(6n**-i)  -  Y,(ank-i  +  bmk-i){tk  -  tk-i) 
fc=i  fc=i  it=i 

(111.16) 

Continuing  with  Step  3,  one  can  take  the  partial  derivatives  of  Equation  III.  16  with 
respect  to  a  and  b  to  get  Equations  III. 17  and  111.18,  respectively. 


da  a       £J 
™ =    -7-  -  Z^mk-i[tk-tk-i)  (111.18) 
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Finally,  maximize  Equations  III.  17  and  III.  18  by  setting  each  of  the  partial  deriva- 
tives equal  to  zero  and  solving  for  a  and  b.  The  resulting  solutions  for  a  and  6, 
denoted  a  and  6,  maximize  the  likelihood  function  (Equation  III.  15)  and  serve  as  the 
COMAN  MLE  estimates  for  the  attrition  rate  coefficients.  See  Equations  III. 19  and  III. 20. 


a  =  ^r — 4 — r^  (IIL19) 

2^k=l  nk-l{tk  —  tk-l) 


J2k=lmk-l{tk  -tk-l) 


2.       Summary  of  the  COMAN  MLE  Method 

While  the  presented  mathematical  description  of  the  COMAN  MLE  method 
for  estimating  attrition  rate  coefficients  is  applicable  only  to  the  homogeneous  rep- 
resentation of  Lanchester's  Square  Law  (F|F)  in  one  of  its  simpler  forms,  the  reader 
can  see  that  the  general  method  is  quite  simplistic.  One  of  the  simplifying  charac- 
teristics of  the  developed  example  is  the  appearance  of  the  attrition  rate  coefficients 
as  linear  terms  only.  It  is  quite  possible  that  the  true  nature  of  the  attrition  rate 
coefficients  is  non-linear  in  actual  combat;  however,  for  the  purposes  of  this  dis- 
cussion, the  author  has  chosen  to  represent  the  mathematics  of  the  method  with 
tractable  forms.  As  a  practical  example,  Appendix  C  details  the  application  of  the 
COMAN  MLE  method  for  attrition  rate  estimation  to  the  high  resolution  output  data 
from  the  AOAC  97-3  Janus  Simulation  described  in  Chapter  II.  The  calculations  are 
presented  for  the  mechanized  battle  between  the  BLUEFOR  M2  (with  TOW  missile) 
and  REDFOR  BMP-2  (with  AT-5  missile).  Note  that  this  computational  excursion 
makes  use  of  Lanchester's  Square  Law  (F|F)  (as  the  mechanized  battle  is  one  char- 
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acterized  by  exchanges  of  aimed  fire)  and  the  mathematical  framework  leading  up  to 
and  including  Equations  III.  19  and  III. 20. 

F.      EXTENSION    OF    ATTRITION   RATE    CONCEPTS 
TO  THE  ATCAL  ALGORITHMS 

By  no  means  should  the  reader  consider  this  treatment  of  attrition  rate  anal- 
ysis to  be  exhaustive.  A  significant  number  of  noteworthy  modifications  have  been 
proposed  and  are  used  within  Lanchester-type  deterministic  models.  For  a  highly 
detailed  inspection  of  these  modifications,  the  interested  reader  should  consult  the 
prolific  works  of  Professor  James  G.  Taylor,  Naval  Postgraduate  School,  and  co- 
advisor  to  this  research.  However,  with  the  general  background  development  of  the 
methods  for  attrition  rate  estimation  complete,  the  reader  is  now  prepared  to  consider 
their  application  to  (or  lack  thereof)  the  ATCAL  algorithms.  The  specific  modifica- 
tions to  the  underlying  Lanchester-type  model  which  form  the  basis  for  the  ATCAL 
methodology  can  be  primarily  attributed  to  Gordon  M.  Clark  and  his  PhD  Thesis, 
"The  Combat  Analysis  Model"  [Ref .  8] ,  while  ongoing  analysis  is  being  conducted  by 
James  G.  Taylor  [Ref.  2].  However,  recall  that  the  specific  attrition  rate  parameter 
estimation  process  used  in  ATCAL  Phase  I  can  be  described  as  an  ad  hoc  process  at 
best. 
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IV.         THE  UNDERLYING 
LANCHESTER-TYPE  MODEL  FOR  ATCAL 

As  the  primary  purpose  of  this  research  is  to  establish  a  test  mechanism  for 
evaluating  the  current  ATCAL  attrition  algorithm  as  it  applies  to  aggregated  combat 
models  (like  CEM  and  TACWAR),  the  interested  reader  must  develop  a  familiarity 
with  ATCAL's  underlying  Lanchester-type  models  and  their  application  via  explicit 
assessment  equations.  As  detailed  documentation  describing  the  theoretical  bases  of 
ATCAL  is  not  available  from  its  original  author,  this  chapter  attempts  to  outline  the 
construct  of  the  ATCAL  attrition  algorithm  based  primarily  upon  a  very  thorough 
treatment  of  the  subject  by  James  G.  Taylor.  [Ref.  2]  It  is  with  an  understanding  of 
the  current  ATCAL  algorithm  that  the  interested  reader  could  conduct  research  for 
improving  its  mathematical  efficiency  or  developing  mathematical  variants  to  fulfill 
other  combat  model  requirements.  The  reader  may  find  it  helpful  to  review  Chap- 
ter III,  Section  C,  which  details  the  basic  framework  of  Lanchester's  Linear  (FT|FT) 
and  Square  (F|F)  Laws  because  of  the  parallelisms  which  extend  to  the  Lanchester- 
type  models  present  in  the  ATCAL  attrition  algorithm. 

A.      THE  HOMOGENEOUS  FORCE  MODEL 

Considering  the  case  of  homogeneous  combat  between  two  opposing  forces,  the 
ATCAL  attrition  algorithm  employs  the  nonlinear  Equations  IV.  1  and  IV. 2  which  can 
be  attributed  to  Gordon  M.  Clark  [Ref.  8]. 


—    =    -a{l  -  px{t))y{t)  with£(0)  =  zo  (IV.l) 

±    =    _^(i_^))xW  withy(0)  =  yo  (IV.2) 


■  (t)  >  0  size  of  force  X  at  time  t 
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y(t)  >  0  size  of  force  Y  at  time  t 

a  >  0  rate  at  which  an  individual  Y  firer  kills  acquired  X  targets 

j3  >  0  rate  at  which  an  individual  X  firer  kills  acquired  Y  targets 

0  <  p       <  1  probability  that  an  individual  Y  firer  has  no 

X  target  available  for  engagement 
0  <  q    '  <  1  probability  that  an  individual  X  firer  has  no 

Y  target  available  for  engagement 

The  reader  may  find  it  more  convenient  to  express  Equations  IV.  1  and  IV. 2 
with  respect  to  target  availabilities  rather  than  non  availabilities.  To  do  so,  let  1  —  (1  — 
^y(t)  _  j  _  p»(*)}  which  is  the  probability  that  an  individual  Y  firer  has  one  or  more 
X  target(s)  available  for  engagement.  And,  analogously,  let  1  —  (1  —  Z?)*W  =  1  —  qy^\ 
With  these  substitutions,  Equations  IV.  1  and  IV. 2  can  be  expressed  as  follows. 

=    _a[i  _  (i  _  A)*W]y(t)  withz(0)=xo  (IV.3) 

at 

^    =    -P[l-(l-ByW]x(t)  w\thy(0)  =  yo  (TV  A) 

Considering  the  elements  of  Equation  IV.3,  the  reader  will  note  that  the  mul- 
tiplication of  the  conditional  kill  rate  (a)  with  the  probability  of  target  availability 
(1  —  (1  —  A)x(*')  produces  the  expected  rate  at  which  an  individual  Y  firer  kills  X 
targets.  Similarly,  the  multiplication  of  (3  and  1  —  (1  —  B)y^  yields  the  expected  rate 
at  which  an  individual  X  firer  kills  Y  targets. 

It  has  been  shown  that  force  levels  can  never  become  zero  or  negative  as  a 
matter  of  theory  (see  [Ref.  9]).  However,  as  the  attrition  calculations  are  necessarily 
executed  on  computer  platforms,  the  potential  for  achieving  "computer  zero"  must 
be  recognized.  So,  in  practice,  the  attrition  models  above  must  be  defined  to  be  valid 
only  when  the  forces  levels,  x(t)  and  y(t)  are  greater  than  zero.  This  refinement  of  the 
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homogeneous  attrition  model  within  the  ATCAL  algorithm  produces  Equations  IV. 5 
and  IV.6. 


dx  -<*[1  -  (1  -  A)*%(t)      for  *(*),</(*)  >0,  m/R. 

It     =     1  ^       5^ 

ai  0  otherwise. 


dy 
dt 


-0[l  -  (1  -  B)«*\x{t)     for  x(t),  y(t)  >  0, 

(IV.6) 

0  otherwise. 


This  positive  force  level  requirement  can  be  easily  incorporated  in  computer  algo- 
rithms by  including  a  conditional  check  of  the  current  force  level  minus  any  calcu- 
lated change  to  the  force  level.  If  the  resulting  force  level  is  less  than  zero,  then  the 
algorithm  should  return  a  value  of  zero.  Otherwise,  the  current  force  level  is  updated 
by  subtracting  the  calculated  attrition. 

1.        Application  of  the  Averaging  Operator 

In  order  to  transform  Equations  IV. 5  and  IV.6  into  explicit  assessment  equa- 
tions as  employed  in  the  ATCAL  algorithm,  it  is  important  to  develop  the  mathe- 
matical properties  of  the  averaging  operator  (over  time)  given  by 

1    r* 
Ave[x]  =  x    =    -       x(s)ds  (IV.  7) 

t  Jo 

where  x(t)  is  a  function  of  time. 

Application  of  this  operator  to  the  Lanchester-type  equations  developed  above  will 
result  in  an  operator  equation  for  the  force  levels,  x(t)  and  y(£),  respectively.  Define 
the  number  of  casualties  inflicted  on  the  X  force  as  Ax  and  the  number  of  casualties 
inflicted  on  the  Y  force  as  Ay. 

Ax    =    x0  —  x(t),         where  x0  =  x(0) 

Ay    =    yo  -  y(t),         where  y0  =  y(0) 

dx 
Now,  consider  applying  the  averaging  operator  to  —  which  results  in 

dt 

.dx,  1    ff  dx 

ds 


dx  1    n  dx 

V6^~dt*  1  Jo   ds 
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=  j[*W-*o] 

[Xp-  X(t)] 
t 

Ax 


t 
Hence,  applying  the  averaging  operator  to  the  homogeneous  Lanchester-type 
Equations  IV. 5  and  IV. 6  representative  of  the  ATCAL  algorithm,  one  obtains  the 
operator  equations,  Equations  IV. 8  and  IV. 9  below. 


J    =    -«[1  -  (1  -  AyW)y(t) 


Ave[^]    =    Ave[-a[l  -  (1  -  AyW]y(i)} 

=    -aAve[[l  -  (1  -  A)x{t)]y(t)] 


and,  similarly, 


t 

Ax    =    atAve[[l  -  (1  -  A)x^}y(t)}  (IV.8) 


Ay    =    (3tAve[[l-{l  -  B)y(t)]x{t)}.  (IV.9) 


The  reader  will  note  that  the  arguments  of  the  averaging  function  appearing  in 
Equations  IV.8  and  IV.9  above  are  highly  nonlinear.  Unfortunately,  these  equations 
cannot  be  further  simplified  to  explicit  expressions  involving  only  average  force  levels 
without  making  some  type  of  approximation.  The  key  approximation  made  within  the 
ATCAL  algorithm  is  that  the  average  of  this  nonlinear  argument  in  x(t)  and  y(t)  is 
equal  to  the  same  functional  relationship  between  x  and  y,  the  results  of  the  averaging 
function  applied  to  x(t)  and  y(t).  Although  this  approximation  is  not  precise,  it  is 
often  applied  to  simplify  averages  of  nonlinear  functional  relationships.  It  is  this 
simplification  that  allows  ATCAL  to  move  forward  computationally.  Equations  IV.8 
and  IV.9  can  now  be  replaced  by  Equations  IV.  10  and  IV.ll. 

Ax    =    at[l  -  (1  -  Af]y  (IV.10) 
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Ay    =    pt[\  -  {I  -  B)*]x  (IV.ll) 

x  =  time  average  of  the  force  level  x(t) 

y  =  time  average  of  the  force  level  y(t) 

Ax  =  number  of  X  casualties  during  time  t 

Ay  —  number  of  Y  casualties  during  time  t 

2.       Applying  the  Results  for  Homogeneous  Forces 

As  written,  Equations  IV. 10  and  IV.ll  have  four  unknown  quantities,  namely, 
the  average  force  levels  (x  and  y)  and  the  numbers  of  casualties  (Ax  and  Ay).  With 
two  equations  and  four  unknowns,  additional  assumptions  are  required  in  order  to 
solve  this  system  of  assessment  equations.  The  second  key  assumption  made  in  the 

ATCAL  algorithm  is  that  average  force  levels  can  be  represented  by  the  expressions 

(  —  Ax)  (—Ay) 

x  =  — -, jr—r  and  y  =  { — .  [Ref.   11    It  so  happens  that  these  expressions 

/nl_Ax  y       /nl_A*)     l  J  vv  F 

hold  true,  respectively,  if  x(t)  and  y(t)  are  exponential  functions  (increasing  or  de- 
creasing) over  time.  As  force  levels  necessarily  decrease  during  battle,  the  negative 
exponential  function  is  assumed.  Combining  this  information  with  Equations  IV.  10 
and  IV.ll,  the  underlying  Lanchester-type  attrition  model  used  to  determine  the 
numbers  of  casualties  during  each  assessment  period  is  represented  in  the  ATCAL 
algorithm  by  the  following  system  of  four  simultaneous  nonlinear  equations  in  the 
four  unknowns,  x,  y,  Ax,  and  Ay. 

Ax    =    erf[l  -  (1  -  Af]y 

Ay    =    0t[l  -  (1  -  Bf}x 

-  (-Ax) 

c  =         (-Ay) 

v       »o  ' 

While  small  perturbations  of  this  notation  are  explicitly  modeled  in  the  ATCAL 
algorithm,  this  system  of  equations  should  be  thought  of  as  the  basic  ATCAL  casualty 
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assessment  system  of  equations.  This  form  better  suits  future  research  in  applying 
numerical  solution  techniques  to  solve  the  system  (recognizing  the  requirement  for 
excluding  negative  force  levels  as  described  in  Equations  IV. 5  and  IV. 6).  The  actual 
method  used  in  ATCAL  to  solve  this  system  of  nonlinear  equations  is  a  heuristic 
form  of  successive  substitutions  which  is  beyond  the  scope  of  this  research.  The 
interested  reader  may  find  the  study  of  alternative  numerical  solution  techniques 
extremely  relevant  to  verifying  current  ATCAL  methodologies  and  proposing  future 
replacement  ones. 

B.   THE  HETEROGENEOUS  FORCE  MODEL 

Considering  the  case  of  heterogeneous  combat  between  two  opposing  forces 
made  up  of  m  distinct  X  force  weapon  systems  and  n  distinct  Y  force  weapon  systems, 
the  ATCAL  attrition  algorithm  employs  the  nonlinear  Clark  Equations  IV.  12  and 
IV.  13  detailed  below.  In  these  equations,  i  =  1  denotes  the  lowest  priority  target 
type  and  i  =  m  denotes  the  highest  priority  target  type  within  the  X  force.  And, 
j  —  I  denotes  the  lowest  priority  target  type  and  j  =  n  denotes  the  highest  priority 
target  type  within  the  Y  force.  Again,  the  reader  may  find  it  helpful  to  compare 
these  equations  with  Equations  III. 8  and  III. 9  of  Chapter  III. 


dxi 
~dt 


71  771 

£  °a(l  "  p5W)    II   PlfvM  for  z  =  1,  2,  •  •  • ,  m 

j=l  k=i+l 


(IV.12) 


J  m  n 

%    =    -5><(1-<#W)    II    «BW*W  fori  =  l,2,..-,n 


dt 


i=l  k=j+l 


(IV.13) 


Xi(t)  >  0  number  of  weapon  systems  of  type  i  in  force  X  at  time  t 

yj(t)  >  0  number  of  weapon  systems  of  type  j  in  force  Y  at  time  t 

£j(0)    =    Xio      for  i  =  1,2,  •  •  •  ,m 
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y;(°)    =    yjo      for  j  =  l,2,---,n 
a,j  >  0  rate  at  which  an  individual  Yj  firer  kills  acquired  X{  targets 

0ji  >  0  rate  at  which  an  individual  X{  firer  kills  acquired  Yj  targets 

0  <  Pij      <  1  probability  that  an  individual  Yj  firer  has  no 

Xi  target  available  for  engagement 
0  <  q*ji      <  1  probability  that  an  individual  Xt  firer  has  no 

Yj  target  available  for  engagement 

Similar  to  the  homogeneous  case,  the  reader  may  find  it  more  convenient  to  express 
Equations  IV. 12  and  IV. 13  in  terms  of  target  availabilities  by  letting  \  —  {l  —  Aij)x%^  = 
1  -  p*j®  and  1  -  (1  -  Bji)y>M  =  1  -  q]l{t).  With  this  substitution,  Equations  IV.12 
and  IV.  13  can  be  expressed  as  follows. 


-^    =    -E^{(l-(l-^)x,(<))   II  (1  -  AkJY^}yi{t)  i  =  l,2,...,m 

at  i=i  k=i+i 

(IV.14) 

j  m  n 

%    =    -E^{(l-(l-^f)^W)    II  C1  -  Bu)»®}xi(t)  j  =  l,2,-..,n 

dt  i=\  k=j+i 

(IV.15) 

Again,  in  order  to  preclude  the  possibility  of  achieving  "machine  zero"  (or 
negative  force  levels)  during  numerical  attrition  calculations,  Equations  IV.14  and 
IV.15  can  be  further  specified  so  that  the  force  levels  X{(t)  and  yj{t)  are  valid  only 
when  they  are  greater  than  zero.  One  final  notational  simplification  renders  the 
heterogeneous  formulation  of  the  underlying  Lanchester-type  model  for  ATCAL  more 
compact.  Consider  the  set  of  indices  corresponding  to  target  types  with  a  higher 
priority  than  a  given  target  type  by  formally  defining 

Uj     =     {k|&  £  integers] 

X  target  type  k  has  a  higher  priority  than  target  type  i  for  a  Yj  firer;  and, 
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Jjj     =     {k\k  €  integers} 

Y  target  type  k  has  a  higher  priority  than  target  type  j  for  a  Xi  firer. 

With  this  simplifying  set  notation,  a  preferred  form  for  the  underlying  Lanchester- 
type  equations  for  ATCAL  appear  below  in  Equations  IV.  16  and  IV.  17.  This  rela- 
tively simple  representation  results  in  less  computational  requirements  when  imple- 
mented as  the  attrition  model  of  a  large-scale  aggregated  combat  model. 


~dt 


-t(0\rT  ^x*(*)' 


-  £?=i  «y{(l  -  $w)Il*ejy  PtjhAt)      ^  *,(*), y3{t)  >  0, 
0  otherwise. 

(IV.16) 

*y        f  -  ££i  M(l  ~  #w)  IW*  «B(tW*)    ^  «,-(*),  y,-W  >  o, 

*"  0  otherwise. 

(IV.17) 

1.       Application  of  the  Averaging  Operator 

Without  loss  of  generality,  the  application  of  the  averaging  operator  (with 
associated  notational  definitions)  to  the  heterogeneous  Lanchester-type  equations  will 
result  in  operator  equations  for  the  force  levels,  £,-(£)  and  yj(t).  In  order  to  maintain 
consistent  notation,  the  author  has  chosen  to  represent  the  equations  using  the  "target 
availability"  notation  of  Equations  IV. 14  and  IV. 15. 


AT.  n                                             m 

i1    =  -£«*{(! -(i-*;)Ii(0)  II  (i  -  4y)x*(t)}y;(0 

dt  j=l                                                               k=i+l 

At-  n                                                              m 

AveQ]     =  ^e[-X>,{(l-(l-^r-W)    nP-^WW)] 

j=l  fc=i+l 

=  -  £«tfiM{(i  -  (l  -  A,r(t))  n  (i  -  4y)*fc(0}y;W] 


*  i=i  fc=*+i 


(IV.18) 
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and,  similarly, 

771  71 

Ay,     =    tJ2^Ave[{(\-(l-BJty^)    I]  (1  -  Bkty^}xt(t)}- 

i=l  k=j+\ 

(IV.19) 

Analogous  to  the  homogeneous  case,  the  key  approximation  made  within  the 
ATCAL  algorithm  is  that  the  average  of  these  nonlinear  arguments  in  x,-(£)  and  yj(t) 
are  equal  to  the  same  functional  relationships  between  x,-  and  y"j,  which  are  the  results 
of  the  averaging  function  applied  to  Xi(t)  and  yj(t).  Applying  this  approximation, 
Equations  IV. 18  and  IV.19  can  be  replaced  by  Equations  IV. 20  and  IV. 21. 


A*,  =  *i)fl%-{(i-(i-^-)*)  na-^rta       (iv.20) 

J=l  k=i+l 

m  n 

Ay,    =    *Eftrf(l-(l -**■)*)   II  (l  -  Bki)Sk}xi  (IV.21) 

t=l  k=j+l 

2.       Applying  the  Results  for  Heterogeneous  Forces 

Similar  to  the  homogeneous  argument  presented  in  the  previous  section,  both 
the  average  force  levels  (x;  and  y^and  the  numbers  of  casualties  (Ax,-  and  Ay,) 
are  unknowns  in  Equations  IV.20  and  IV.21.  Therefore,  this  system  is  made  up  of 
(m  +  n)  nonlinear  equations  in  2  •  (m  +  n)  unknowns.  In  order  to  increase  the  number 

of  equations  to  equal  the  number  of  unknowns  ,  ATCAL  again  makes  use  of  an 

(-Axi) 
assumption  for  the  average  force  levels  resulting  in  the  expressions  X{  =  — ^— ■ 

V  X.n    / 


(-Ay,) 


a^.O 


and  yj  =  -r- — .    [Ref.    11]  Combining  this  information  with  Equations  IV.20 

ln{\ u) 

and  IV.21,  the  underlying  Lanchester-type  attrition  model  used  to  determine  the 
numbers  of  casualties  during  each  assessment  period  is  represented  in  the  ATCAL 
algorithm  as  the  following  system  of  2  •  (m  +  n)  simultaneous  nonlinear  equations  in 
the  2  •  (m  +  n)  unknowns,  £;,  yj,  Ax;,  and  Ay.,-. 
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Ax,  =  *2^o^{(i-(i-i^0'0r&H.i(i-^)-*}yi 

Ay,    =    t  Si  fti{(l  -  (1  -  ^)*)IlLi+i(l  ~  BK)fc}*.- 

(-Ax,) 
/n(l-^i) 

While  small  perturbations  of  this  notation  are  explicitly  modeled  in  the  ATCAL 
algorithm,  this  system  of  equations  should  be  thought  of  as  the  basic  ATCAL  caualty 
assessment  equations.  This  form  better  suits  future  research  in  applying  numerical 
solution  techniques  to  solve  the  system  (recognizing  the  requirement  for  excluding 
negative  force  levels  as  described  in  Equations  IV.  16  and  IV.  17).  The  actual  method 
used  in  ATCAL  to  solve  this  system  of  nonlinear  equations  is  a  heuristic  form  of 
successive  substitutions  which  is  beyond  the  scope  of  this  research.  The  interested 
reader  may  find  the  study  of  alternative  numerical  solution  techniques  extremely 
relevant  to  verifying  current  ATCAL  methodologies  and  proposing  future  replacement 
ones. 

C.      SUMMARY  OF  THE  LANCHESTER-TYPE  ASSESS- 
MENT EQUATIONS  IN  ATCAL 

The  homogeneous  and  heterogeneous  systems  of  nonlinear  equations  found 
at  the  ends  of  Sections  A  and  B,  respectively,  comprehensively  describe  the  under- 
lying Lanchester-type  models  and  their  basic  assessment  equations  for  the  ATCAL 
attrition  methodology.  "The  most  significant  implication  of  these  results,  however, 
is  that  a  much  more  convenient  way  of  generating  numerical  results  for  ATCAL 
has  immediately  been  established:  namely,  numerical  integration  of  the  underlying 
Lanchester-type  differential  equations."  [Ref.  2]  Moreover,  with  these  mathematical 
constructs,  the  interested  reader  is  now  prepared  to  further  investigate  the  soundness 
of  the  heuristic  successive  substitution  method  of  ATCAL  and  compare  it  to  more 
robust  numerical  solution  techniques. 
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V.         CONCLUSION 

A.  CONCLUSIONS 

This  research  has  developed  a  practical  test  mechanism  for  studying  the  attri- 
tion methodologies  of  aggregated  combat  models.  The  Janus  high  resolution  simula- 
tion software,  capable  of  modeling  brigade  sized  (and  below)  operations,  is  completely 
adequate  for  producing  the  required  high  resolution  output  data  required  for  the  anal- 
ysis of  aggregated  attrition  methodologies.  Janus  is  both  widely  available  within  the 
military  community  (to  include  the  Naval  Postgraduate  School)  and  is  easily  learned 
through  the  use  of  a  tutorial.  See  [Ref.  5].  This  user  friendly  computer  simulation 
provides  the  researcher  with  sufficient  detail  and  a  mechanism  for  capturing  detailed 
interactions  on  the  modeled  battlefield.  The  ability  to  capture  the  highly  detailed 
interactions  between  weapon  systems  is  what  qualifies  Janus  as  a  valuable  tool  for 
conducting  analysis  of  the  attrition  algorithms  inherent  to  ATCAL  (and  other  ag- 
gregated attrition  methodologies).  The  reader  should  follow  the  steps  outlined  in 
Section  C  to  initiate  such  a  study. 

B.  RECOMMENDATIONS 

First,  having  established  the  vehicle  (Janus)  for  generating  the  data  required 
for  developing  attrition  (kill)  rate  parameters,  the  author  has  opened  the  door  for 
continued  research  in  the  area  of  attrition  methodologies.  As  this  research  developed 
and  discussed  one  such  combat  scenario,  the  continued  study  of  attrition  calculations 
would  require  additional  Janus  scenarios  with  varying  parameters.  The  creation  and 
execution  of  such  scenarios  could  easily  be  incorporated  into  current  high  resolution 
modeling  courses  taught  at  the  Naval  Postgraduate  School. 

Second,  as  discussed  in  Chapter  IV,  the  nonlinear  attrition  assessment  equa- 
tions resulting  from  the  Lanchester-type  equations  used  in  the  ATCAL  algorithm  are 
currently  "solved"  using  a  heuristic  method  of  successive  substitutions.    Significant 
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follow-on  research  is  applicable  in  the  employment  of  a  more  robust  numerical  so- 
lution technique  for  these  systems  of  equations.  This  area  of  continued  study  does 
not  necessarily  depend  upon  additional  Janus  scenarios;  however,  additional  scenarios 
would  provide  a  data  source  for  technique  verification. 

Third,  the  computer  code  which  contains  the  ATCAL  Phase  I  and  II  algo- 
rithms exists  in  several  mediums.  A  significant  endeavor  integral  to  making  ATCAL 
(and  any  proposed  modifications)  more  efficient  includes  a  thorough  analysis  of  the 
current  computer  code  which  exists  in  at  least  two  languages,  FORTRAN  and  C. 
The  interested  reader  may  find  that  the  algorithms,  which  make  up  ATCAL  would 
be  more  efficient  in  one  or  more  specific  computer  languages.  For  example,  a  proposed 
numerical  solution  technique  may  be  more  conducive  to  one  computer  language  than 
another. 

Finally,  an  automated  bridge  is  required  between  the  Janus  high  resolution 
output  data  and  the  computer  code  used  to  perform  the  aggregated  attrition  calcula- 
tions (like  ATCAL).  While  this  capability  may  not  be  directly  applicable  to  the  actual 
execution  of  aggregated  models  as  they  currently  exist,  it  is  integral  to  the  efficient 
study  of  the  algorithms  themselves.  It  is  in  the  research  arena  that  this  automated 
capability  for  linking  Janus  output  to  attrition  algorithms  would  be  helpful.  JETS 
may  be  an  important  tool  for  continued  study  of  this  recommendation. 

As  combat  models  are  an  integral  part  of  the  military  analysis  and  training 
communities,  continued  analysis  of  attrition  methodologies  (specifically  ATCAL)  in 
the  field  of  aggregated  combat  modeling  will  greatly  enhance  model  believability  and 
overall  usefulness.  As  computer  technologies  continue  to  outpace  attrition  algorithm 
development,  verification,  and  validation,  a  continued  focus  towards  efficiency  and 
accuracy  will  ensure  the  applicability  of  current  and  future  attrition  methodologies. 
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C.      INITIATION    OF    ATTRITION    ANALYSIS    USING 
JANUS 

The  key  to  initiating  a  detailed  study  of  attrition  and  its  inherent  algorithms 
lies  within  the  generation  of  high  resolution  combat  simulation  data.  As  such,  the 
researcher  must  carefully  consider  the  process  of  creating  high  resolution  scenarios. 
The  interested  researcher  should  consider  the  following  steps  which  outline  the  model 
formulation  and  data  collection  process. 


Step  1 

Define  the  tactical  scenario  (forces  composition,  terrain  type, 
environmental  conditions,  and  concept  of  the  operation) 
for  which  a  model  is  desired. 

Step  2 

Construct  the  Janus  simulation  scenario  and  validate  its 
representation  of  the  desired  interaction  between  modeled  forces. 

Step  3 

Determine  an  appropriate  number  (n)  of  iterations  for  which  to 
execute  the  Janus  simulation. 

Step  4 

Independently  execute  the  Janus  simulation  n  times  to  develop 
a  robust  sample  of  simulation  output  data  from  which  statistically 
significant  killer/victim  information  can  be  determined. 

Step  5 

Analyze  the  consolidated  high  resolution  output  data  and  estimate 
the  attrition  (kill)  rates  between  opposing  weapons  systems. 

Step  6 

Incorporate  the  derived  attrition  (kill)  rates  in  the  analysis  of 
attrition  algorithms  inherent  to  ATCAL  (and  other  aggregated 
attrition  methodologies). 

Step  7 

Utilize  an  appropriate  Lanchester-type  deterministic  model  to 
replay  the  simulated  battle. 

Table  XI.  Steps  for  Attrition  Analysis  and  Deterministic  Battle  Replay  Using  Janus 
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APPENDIX  A.  JANUS  COMBATANTS 
LISTING  AND  RUN-TIME  BATTLEFIELD 

PICTURES 

This  Appendix  provides  the  reader  with  the  detailed  description  of  the  combat 
systems  making  up  the  REDFOR  and  BLUEFOR  modeled  in  the  AOAC  97-3  Janus 
Simulation  described  in  Chapter  II.  Additionally,  pictures  are  included  to  depict  the 
battlefield  details  and  force  allocations  prior  to  the  start  of  hostilities  and  at  30  minute 
increments  throughout  the  simulation. 

T-80  Main  Battle  Tank  15 

BMP-2  Armored  Personnel  Carrier  55 

BRDM-A  Armored  Personnel  Carrier  10 

BRDM-2  Armored  Personnel  Carrier  7 
GMZ-MI  Armored  Mine  Warfare  System     2 

MTK.LI  Armored  Engineer  System  1 

2S12  M  122mm  Self-Propelled  Artillery  16 

IF  TRU  Wheeled  Truck  26 

DF  Tru  Wheeled  Truck  16 

2S1  155mm  Self-Propelled  Artillery  30 

BM-21  Rocket  Launcher  8 

SNAR-2  Radar  2 

ZOOPAR  Radar  2 

2S5  152mm  Self-Propelled  Artillery  36 

ADA  TM  Air  Defense  Artillery  Team  5 

ZSU-23  Air  Defense  Artillery  3 

HIND  Attack  Helicopter  6 

POL  AI  Aviation  Fuel  Depot  2 


Total  REDFOR  Combatants  242 

Table  XII.  REDFOR  Organization  for  AOAC  97-3  Janus  Simulation 
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M1A1  Main  Battle  Tank  82 

M1A1P  Main  Battle  Tank  w/  Mine  Plow  24 

Ml  AIR  Main  Battle  Tank  w/  Mine  Roller  8 

M2  Armed  Armored  Personnel  Carrier  113 
BSFV  Armed  Improved  Armored  Personnel  Carrier     6 

Ml  13  Armored  Personnel  Carrier  21 

FIST-V  Artillery  Command  Vehicle  4 

CEV  Combat  Engineer  Vehicle  4 

BRIDGE  Mobile  Bridge  Carrier  5 

ADA  TE  Air  Defense  Artillery  Team  4 

ACE  Armored  Combat  Earthmover  13 

AVLM  Armored  Vehicle  Launched  Bridge  9 

AVLM-W  Armored  Vehicle  Launched  Bridge  6 

AVLB  W  Armored  Vehicle  Launched  Bridge  2 

HMV.50  Wheeled  Vehicle  47 

HMV/M1  Wheeled  Vehicle  10 

HMV/SM  Wheeled  Vehicle  6 

HMV/CO  Wheeled  Vehicle  8 

IF  TRU  Wheeled  Truck  25 

DF  Tru  Wheeled  Truck  4 

M106A1  120mm  Mortar  Vehicles  22 

M577  Armored  Personnel  Carrier  (Command)  6 

GSR/Ml  Ground  Surveillance  Radar  1 

SCT  TE  Scout  Team  34 

SAPPER  Combat  Engineer  Team  90 

SAW  TE  Machine  Gun  Team  144 

VOLCANO  Mine  Dispensing  System  4 

Ml 09 A3  155mm  Self-Propelled  Artillery  24 

M109A6  155mm  Self-Propelled  Artillery  18 

HMMQ36  Artillery  Radar  System  1 

ADA  TE  Air  Defense  Artillery  Team  4 

SMK113  Armored  Vehicle  Smoke  Generator  3 

A- 10  Fixed  Wing  Close  Air  Support  12 

Total  BLUEFOR  Combatants  764 

Table  XIII.  BLUEFOR  Organization  for  AOAC  97-3  Janus  Simulation 
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Figure  1.  BLUEFOR  Initial  Disposition  for  AOAC  97-3  Janus  Simulation 
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Figure  2.  REDFOR  Initial  Disposition  for  AOAC  97-3  Janus  Simulation 
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Figure  3.  Battlefield  Picture  at  Start  (0  Minutes)  for  AOAC  97-3  Janus  Simulate 


;ion 
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Figure  4.  Battlefield  Picture  at  30  Minutes  for  AOAC  97-3  Janus  Simulation 
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Figure  5.  Battlefield  Picture  at  60  Minutes  for  AOAC  97-3  Janus  Simulation 
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Figure  6.  Battlefield  Picture  at  90  Minutes  for  AOAC  97-3  Janus  Simulati 


ion 
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Figure  7.  Battlefield  Picture  at  120  Minutes  for  AOAC  97-3  Janus  Simulation 
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Figure  8.  Battlefield  Picture  at  150  Minutes  for  AOAC  97-3  Janus  Simulation 
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Figure  9.  Battlefield  Picture  at  180  Minutes  for  AOAC  97-3  Janus  Simulation 
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Figure  10.    Battlefield  Picture  at  Completion  (211  Minutes)  for  AOAC  97-3  Janus 
Simulation 
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APPENDIX  B.  APPLICATION  OF  THE 

ANDERSON  FR  METHOD  TO  A 

HYPOTHETICAL  HETEROGENEOUS 

BATTLE  REPRESENTATIVE  OF  THE 

AOAC  97-3  JANUS  SIMULATION 

This  appendix  provides  the  reader  with  an  example  application  of  the  Ander- 
son FR  Method  with  subordinate  APP  process  for  attrition  rate  estimation  using 
hypothetical  kill  rate  data  representative  of  the  simulation  output  of  the  AOAC  97-3 
Janus  Simulation  described  in  Chapter  II.  Specifically,  the  author  has  selectively 
chosen  to  calculate  the  BLUEFOR  and  REDFOR  Firepower  Scores  (FPS)  and  the 
resultant  Firepower  Index  (FPI)  using  two  unique  BLUEFOR  weapon  systems  and 
three  unique  REDFOR  weapon  systems.  The  heterogeneous  aimed  fire  battle  between 
these  opposing  weapons  systems  provides  the  reader  with  a  numerically  obvious  closed 
form  representation  for  both  forces'  FPS. 

1.       FORCES  COMPOSITION  AND  NOTATIONAL  TRANS 
LATION 

Let  X  represent  the  BLUEFOR  composed  of  two  weapon  system  types,  namely, 
114  Ml  Main  Battle  Tanks  and  113  M2  Bradley  Fighting  Vehicles.  Let  Y  represent 
the  REDFOR  composed  of  three  weapon  system  types,  namely,  15  T-80  Main  Battle 
Tanks,  55  BMP-2  Armored  Personnel  Carriers,  and  17  BRDM  Armored  Personnel 
Carriers.  Recall  from  Chapter  III  that  the  goal  of  the  Anderson  FR  Method  is  the 
representation  of  the  heterogeneous  force  as  a  single  scalar  measure  of  combat  power. 
Translating  Equation  III.l  from  Chapter  III,  we  have  the  mathematical  expressions 
detailed  below  by  Equations  B.l  and  B.2. 
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FPIX  =  J^SfXi  (BLUEFOR  Firepower  Score)  (B.l) 

i  =  integer  index  where  1  =M1  and  2  =M2 

Si  =  value  (score)  of  one  Ml 

S*  =  value  (score)  of  one  M2 

Xi  =  114  Mis 

X2  =  113  M2s 


3 

FPIy    =    52S?Yj  (REDFOR  Firepower  Score)  (B.2) 

i=i 
j    =    integer  index  where  1  =T-80,  2  =BMP-2,  and  3  =BRDM 

S^  =  value  (score)  of  one  T-80 

52  =  value  (score)  of  one  BMP-2 

53  =  value  (score)  of  one  BRDM 
Yi  =  15  T-80s 

y2    =    55  BMP-2s 
Y3    =    17BRDMs 

2.       DETERMINATION    OF   WEAPON    SYSTEM    VAL- 
UES 

Recall  that  the  fundamental  methods  for  determining  the  values  of  Sf  and 
Sj  include  the  analysts'  perception  of  combat  value,  historical  combat  performance, 
measures  of  weapon  firepower  and  support  platform  characteristics,  and,  measures  of 
what  a  weapon  can  kill.  The  last  method  attempts  to  compute  the  value  (score)  for 
a  weapon  as  being  proportional  to  the  rate  at  which  it  destroys  the  value  of  opposing 
enemy  systems.  This  method  for  determining  weapons  system  scores,  called  the 
AntiPotential-Potential  (APP)  method,  is  the  primary  focus  of  this  appendix. 
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The  formulation  of  the  APP  problem  relies  upon  the  "rate"  at  which  systems 
destroy  opposing  systems  in  battle.  These  assumed  or  derived  non-negative  scalar  at- 
trition rates  (one  for  each  weapon  system  firing  at  every  opposing  weapon  system  for 
both  forces)  account  for  the  highly  detailed  engagement  process  described  in  Chap- 
ter I.  While  the  actual  kill  rates  for  each  pairwise  combination  of  weapon  systems 
resulting  from  the  Janus  simulation  described  in  Chapter  II  could  be  derived  through 
analyzing  the  killer/victim  scoreboard,  the  author  has  chosen  to  use  hypothetical  kill 
rate  values  for  each  pair  of  interacting  weapon  systems  thus  defining  kill  rate  ma- 
trices, Kij  and  Lji  below.  This  simplification  assists  the  reader  with  focussing  on 
the  task  at  hand,  understanding  the  Anderson  FR  methodology  and  the  subordinate 
APP  process.  An  example  derivation  of  attrition  rates  from  high  resolution  output 
is  discussed  in  detail  in  Appendix  C. 

.4    .5 


KX3  = 


A    .5    .5 
.2    .3    .4 


and 


Lji  — 


.3    .4 
.2    .3 


The  matrix  Kij  represents  the  rate  at  which  one  type  i  weapon  system  (1=M1  and 
2=M2)  kills  type  j  weapon  systems  (l=T-80,  2=BMP-2,  and  3=BRDM).  Analo- 
gously, Lji  represents  the  rate  at  which  one  type  j  weapon  system  kills  type  i  weapon 
systems.  As  an  illustration,  the  reader  will  note  that  the  BLUEFOR  (X  force)  Ml 
Main  Battle  Tank  kills  the  REDFOR  (Y  force)  BMP-2  Armored  Personnel  Carrier 
at  the  rate  of  K12  —  -5. 

The  remaining  unknowns  of  Equations  B.l  and  B.2  are  the  weapon  systems 
scores,  S*  and  Sj .  Deriving  these  individual  scores  equates  to  solving  Equations  B.3 
and  B.4  (below)  for  the  unknowns  Cx,Cy,S* ,  and  Sj . 


CxSf    =    'EKijSj  for  1  =  1,2 

2=1 

CySj    =    j^LjiSf  forj  =  1,2,3 


(B.3) 

(B.4) 


i=\ 
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Cx     —    proportionality  constant 
Cy    =    proportionality  constant 

So,  for  example,  the  value  of  the  BLUEFOR  Ml  Main  Battle  Tank  (S*)  is  propor- 
tional to  (Cx)  the  sum  total  value  of  every  REDFOR  system  that  it  kills.  The  reader 
will  note  that  the  first  equation  above  represents  two  equations  (one  each  for  i  =  1,2) 
and  the  second  equation  represents  three  equations  (one  each  for  j  =  1,2,3)  for  a 
total  of  five  equations.  Moreover,  this  system  of  five  equations  has  seven  unknown 
variables,  Cx,  Cy,  S* ^S* -,  Sf ,  S% ,  and  S% .  A  system  of  linear  equations  with  more 
unknowns  than  equations  defining  them  is  referred  to  as  underdetermined. 

a.        Developing  the  Representative  Eigenvalue  Prob- 
lem 

To  solve  the  value  equations  defined  by  Equations  B.3  and  B.4,  consider  the 

following  substitutions. 

CySj    =    £>;,S?  forj  =  1,2,3 

^SJ    =    TriJl^Sf)  for  3  =  1,2,3 

CY    l=l 

Substituting  the  expression  for  Sj  above  into  Equation  B.3  yields 

CxSf    =    J2K»[^-(t^Sf)]  fori  =  1,2 

j=l  ky    i=1 

2       3 

^CxCyS?    =    EEWi*5?)  fori  =  1,2,  (B.5) 

and,  applying  the  same  construct  to  Equation  B.4 


3        2 

CxCySj    =    EEfe^?)  for  j  =  1,2,3.  (B.6) 

3=1  i-\ 
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For  clarity,  Equations  B.5  and  B.6  can  be  expressed  more  compactly  us- 
ing vector  and  matrix  notation  as  depicted  by  Equations  B.7  and  B.8.    Note  that 
Iw=i  T%=i{KijLji)  =  {KL)  is  a  square  matrix  of  dimension  (2x2)  and  Yl%i  T%=i(LjiKij) 
(LK)  is  a  square  matrix  of  dimension  (3x3). 

f  .2600    .3500    .4000 


KL  = 


.4100  .5500 
.2500  .3400 


and 


LK  = 


.2000  .2700  .3100 
.1400  .1900  .2200 


CXCYSX    =    (KL)SX 
CxCyS*    =    (LK)^ 


(B.7) 
(B.8) 


By  letting  A  =  CxCy  [Ref.  9],  the  reader  will  note  that  Equations  B.7  and  B.8 
can  be  written  as  depicted  below  and  represent  a  pair  of  eigenvalue  problems  which 
will  share  common  eigenvalues. 


XSX    =    [KL)SX 
XS*    =    (LK)§* 


(B.9) 
(B.10) 


The  reader  should  note  that  the  matrix  products  (KL)  and  {LK)  are  square 
matrices  which  have  only  positive  entries.  Consequently,  the  Perron- Fro benius  The- 
orem states  that  the  largest  eigenvalue  of  a  matrix  with  only  positive  entries  and 
the  entries  of  the  eigenvector  associated  with  the  largest  eigenvalue  are  positive  real 
numbers.  See  [Ref.  12]. 

b.       Solving  the  Representative  Eigenvalue  Problem 

The  eigenvalue  problems  represented  by  Equations  B.9  and  B.10  can  be  solved 
by  finding  the  roots  of  the  respective  characteristic  equations, 

det(KL-XI)    =    0 
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and, 


det( 


.4100  -  A   .5500 
.2500   .3400  -  A 


)    =    o 


(B.ll) 


det{LK-\I)  0 

.2600  -  A        .3500  .4000 

det{        .2000        .2700 -A        .3100        )  0 

.1400  .1900        .2200  -  A 

where  det  is  defined  as  the  determinant  and  /  is  the  identity  matrix  of  the  appropriate 
dimension.  These  characteristic  equations  are  simplified  to  the  following  equations 
by  evaluating  the  determinant  and  collecting  like  terms.  If  required,  the  reader  can 
review  the  methods  for  evaluating  determinants  of  (2x2)  and  (3x3)  matrices  in  [Ref. 
12]. 


det(KL  -  XI)    =    (.4100  -  A)(.3400  -  A)  -  (.2500)(.5500) 
=    A2  -  .7500A  +  .0019 

and 

det{LK-XI)    =    (.2600-A)[(.2700-A)(.2200-A)-(.1900)(.3100)] 
-.3500[(.2000)(.2200  -  A)  -  (.1400)(.3100)] 
+.4000[(.2000)(.1900)  -  (.1400)(.2700  -  A)] 

=    -A3  +  .7500A2  -  .0019A 

=    -A(A2  -  .7500A  +  .0019) 

Setting  both  determinants  equal  to  zero  and  solving  for  the  values  of  A  which 
fulfill  these  equalities  (using  the  quadratic  formula  [Ref.   12]),  we  find  that  the  first 
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characteristic  equation  yields  roots  Xt  =  .7475  and  A2  =  .0025,  and,  the  second 
equation  yields  Xi  =  .7475,  A2  =  .0025,  and  A3  =  0.  The  reader  is  now  able  to  confirm 
the  existence  of  the  largest  eigenvalue,  A*  =  .7475,  which  is  both  real  and  positive 
as  ensured  by  the  Perron- Fro benius  Theroem.  Substituting  this  value,  A*  =  .7475, 
into  Equations  B.9  and  B.10  presents  the  following  eigenvector  problems  for  which 
we  seek  the  solutions. 


{KL-X*I)SX    =    0 
(LA'-A*/)5y    =    0 

Consider  Equation  B.12  in  isolation. 


(B.12) 
(B.13) 


(KL-X*I)SX    =    0 


Sx 


S? 


-.3375      .5500 
.2500      -.4075 

As  there  are  infinitely  many  eigenvectors  which  solve  this  equality,  the  author  will 
arbitrarily  fix  Sx  =  1  to  produce  a  unique  value  for  Sx  (and  hence,  a  unique  eigen- 
vector). This  amounts  to  normalizing  the  value  of  the  BLUEFOR  Ml  Main  Battle 
Tank  equal  to  one,  which  will  force  all  of  the  weapon  systems  values  to  be  scaled 
relative  to  the  Ml.  The  reader  should  be  aware  that  this  scaling  method  is  only  one 
from  a  list  of  many  methods.  Further  information  can  be  found  in  [Ref.  9].  Solving 
for  the  second  unknown  (the  value  of  the  BLUEFOR  M2  Bradley  Fighting  Vehicle), 

the  reader  can  confirm  that  Sx  =  .6134.  Therefore,  Sx  = 

However,  the  reader  should  be  aware  of  the  fact  that  this  solution  vector  is 
not  unique.  Any  scalar  multiple  of  this  value  vector  will  also  satisfy  Equation  B.12. 
The  proportionality  constants,  Cx  and  Cy ',  appearing  in  Equations  B.7  and  B.8  will 
simply  adjust  to  absorb  the  scalar  multiple  changes  to  the  eigenvector  of  weapon 
system  values  -  a  condition  for  which  the  chosen  scaling  method  shortly  corrects. 


.6134 
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Now  consider  the  effectiveness  of  the  BLUEFOR  Ml  Main  Battle  Tanks  as 
they  destroy  the  "value"  of  the  REDFOR  weapon  systems.   This  can  be  evaluated 

3 

mathematically  by  the  expression  ^^KijSj.    The  ID  AG  AM  scaling  method  [Ref. 

3 

9]  seeks  to  scale  the  values  of  Sj  such  that  ^  KijSj  =  vA*,  where  A*  =  .7475 
(the  largest  eigenvalue).  Recall  that  in  the  development  of  the  eigenvalue  problems 
represented  by  Equations  B.9  and  B.10  that  the  author  let  A  =  CxCy.   Combining 

3 

Equation  B.3  with  the  known  information  S*  =  1  and  ^^KijSj  =  ^^  ^e  rea(^er 
will  notice  the  following  result. 


3=1 

cx  =  >A* 


Cx    =    yCx  Cy 

(Cx)       =     CxCy 
Cx     =     Cy  =  vA* 


Substituting  the  kill  rate  information  from  matrix  K  into  the  equation  ^  Kij  Sj   = 

3=1 

vA*  results  in  the  equation  .45^  +  .55^  +  .5S%  =  .8646.  This  particular  scaling 
method  now  calls  for  the  replacement  the  first  equation  of  the  linear  system  repre- 
sented by  Equation  B.13  with  this  new  equation.  This  action  results  in  the  following 
matrix  equation. 


[LK-yi)^  0 

.4000  .5000  .5000 
.2000  -.4775  .3100 
.1400      .1900      -.5275 


\sr' 

.8646 

si 

= 

0 

si 

0 
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Several  methods  exist  for  solving  this  system  of  three  linear  equations  in  three  un- 
knowns (many  of  which  are  detailed  in  [Ref.    12]  and  [Ref.    13]).   When  solved,  the 

.8176 
resulting  values  for  the  REDFOR  weapon  systems  are  S    =      .6309 

.4442 


3.       CALCULATING    THE    RESULTANT    FIREPOWER 
INDICES 

Armed  with  the  values  for  each  of  the  BLUEFOR  and  REDFOR  distinct 
weapon  systems,  the  reader  can  easily  calculate  the  single-valued  scalar  descriptors 
FPIX  and  FPIy. 


FPIX    =    52S?Xi 

=    (1)(114)  +  (.6134)(113) 
=    183.3142 


=    (.8176)(15)  +  (.6309)(55)  +  (.4442)(17) 
=    54.5151 

As  described  in  Chapter  III,  the  Force  Ratio  (FR^-y)  can  be  calculated  as  FR^y  = 

=  =  3.3626.    In  a  simplified  sense,  this  resultant  FRj^y  would  de- 

rrly  54.5151 

termine  the  outcome  of  a  battle  with  respect  to  parameters  set  within  an  aggregated 
model.  For  a  more  detailed  treatment  of  the  theoretical  foundations  of  the  Anderson 
(FR)  Method  and  a  critique  of  its  use  in  force-ratio-based  attrition  models,  see  [Ref. 
9]. 

The  reader  may  have  gathered  that  the  complexity  associated  with  hetero- 
geneous battles  made  up  of  many  distinct  weapon  systems  would  require  the  intro- 
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duction  of  an  automated  capability  to  solve  the  representative  eigenvalue  problems. 
For  example,  thoroughly  accounting  for  the  33  BLUEFOR  and  18  REDFOR  distinct 
weapon  systems  modeled  in  the  Janus  simulation  described  in  Chapter  II  would  result 
in  a  linear  system  of  51  equations  in  53  unknowns  (counting  Cx  and  Cy).  The  reader 
may  also  surmise  that  the  kill  rate  matrices  associated  with  this  complex  simulation 
may  be  sparse;  in  other  words,  many  of  the  entries  of  the  kill  rate  matrices  could  be 
equivalent  to  zero  based  upon  weapon  systems  which  do  not  directly  interact  on  the 
battlefield.  This  fact  introduces  another  consideration  to  the  Anderson  FR  method 
and  the  subordinate  APP  process  —  the  possibility  for  ill-conditioned  matrices  (see 
[Ref.  12]).  At  this  point,  the  reader  should  have  a  basic  understanding  of  the  An- 
derson (FR)  based  attrition  model  and  its  application  to  a  simplified  heterogeneous 
battle. 
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APPENDIX  C.  APPLICATION  OF  THE 

COMAN  MAXIMUM  LIKELIHOOD 

ESTIMATION  METHOD  TO  A 

HOMOGENEOUS  BATTLE  WITHIN  THE 

AOAC  97-3  JANUS  SIMULATION 

This  appendix  provides  the  reader  with  an  example  application  of  the  COMAN 
MLE  method  for  attrition  rate  estimation  using  portions  of  the  simulation  output 
data  from  the  AOAC  97-3  Janus  Simulation  described  in  Chapter  II.  Specifically,  the 
author  has  chosen  to  apply  this  method  to  the  time  series  of  attrition  between  two 
opposing  mechanized  infantry  weapon  systems.  The  aimed  fire  battle  between  the 
BLUEFOR  M2  Bradley  Fighting  Vehicles  and  the  REDFOR  BMP-2  Armored  Per- 
sonnel Carriers  is  examined  as  a  homogeneous  exchange  of  direct  fires.  The  primary 
weapon  systems  engaged  in  this  battle  are  the  TOW  anti-tank  missile  (BLUEFOR) 
and  the  AT-5  anti-tank  missile  (REDFOR).  As  such  the  derived  attrition  rate  pa- 
rameters will  be  applicable  to  the  Lanchester  Square  Law  (F|F)  for  aimed  fire  war- 
fare. Finally,  the  author  develops  a  practical  application  of  the  differential  equations 
(Equations  C.l  and  C.2  detailed  below)  to  validate  the  derived  attrition  rate  esti- 
mates. This  replay  is  analogous  to  the  replay  executed  within  the  ATCAL  Phase  II 
algorithm.  (However,  the  reader  should  recall  that  ATCAL  Phase  I  uses  an  ad  hoc 
method  for  attrition  rate  estimation  rather  than  the  COMAN  MLE  method.)  The 
data  for  this  analysis  was  exported  from  the  Janus  simulation  output  file  and  manipu- 
lated in  a  commercial  spreadsheet  software  package  (not  included  in  this  Appendix). 
The  data  appears  in  Table  XIV  on  page  104  of  this  Appendix.  However,  the  in- 
terested reader  can  obtain  copies  of  the  associated  spreadsheet  files  from  Professor 
Bard  K.  Mansager,  Department  of  Mathematics,  Naval  Postgraduate  School. 
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1.       NOTATION 

Let  X  represent  the  BLUEFOR  M2  systems  (with  TOW  missiles)  and  Y 
the  REDFOR  BMP-2  systems  (with  AT-5  missiles).  The  relative  Lanchester  Square 
Law  equations,  modified  with  the  COMAN  MLE  attrition  rate  estimates,  model  the 
change  in  force  size  with  respect  to  time.  These  are  Equations  C.l  and  C.2. 


dx 

—     =    -ay(t)  with  x(0)  =  x0  (C.l) 

at 

^     =    ~bx(t)  withy(0)  =  yo  (C.2) 


x(t)  =  size  of  force  X  at  time  t 

y(t)  =  size  of  force  Y  at  time  t 

a  =  derived  attrition  rate  at  which  BMP-2  firers  kill  M2  targets 

b  =  derived  attrition  rate  at  which  M2  firers  kill  BMP-2  targets 

Recall  Equations  III. 21  and  III. 22  on  page  44  of  Chapter  III,  as  shown  below. 


a    — 


J2k=i  nk-i{tk  —  tk-i) 


rY 
b     = 


Ejfc=l  mk-l(tk  —  tk-l) 

2.       DEVELOPMENT  OF  THE  ATTRITION  RATE  ESTI- 
MATES 

Applying  Equations  III. 21  and  III. 22  to  the  Janus  output  data  found  in  Ta- 
ble XIV,  we  have  the  following: 
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a    = 


Hk=l  nk-l(tk  —  tk-l) 

Total  No.  of  M2  Casualties  (K) 
XZ^_j(No.  of  BMP-2  Surviving  at  tk  —  1)  •  (Time  between  kth  and  k  —  lst  Casualty) 

17 
469493  seconds 
=    0.00217256  M2  kills  per  BMP-2  firer  per  minute;       and, 


6    = 


CT 


EjELi  rnk-\{tk  -  tk-i) 

Total  No.  of  BMP-2  Casualties  (K) 
Hj[Li(No.  of  M2  Surviving  at  tk  —  1)  •  (Time  between  kth  and  k  —  lst  Casualty) 
12 
~     997883  seconds 
=    0.00072153  BMP-2  kills  per  M2  firer  per  minute. 

Intuitively,  from  the  standpoint  of  classic  military  operations,  the  resulting 
attrition  rate  estimates  are  logical  with  respect  to  the  modeled  combat  scenario.  The 
REDFOR  is  executing  a  layered  defense  with  fortified  positions  and  established  fields 
of  fire,  while  the  BLUEFOR  is  on  the  attack  characterized  by  weapon  system  exposure 
and  less  distinct  fields  of  fire.  Therefore,  as  one  would  expect,  the  REDFOR  BMP-2 
kills  the  BLUEFOR  M2  at  a  higher  rate  than  vice  versa  (all  else  being  equal).  This 
provides  the  analyst  with  a  common  sense  oriented  way  of  justifying  his  or  her  inter- 
mediate results. 

3.       DERIVATION  OF  FORCE  LEVEL  EQUATIONS  FROM 
LANCHESTER'S  SQUARE  LAW  (F|F)  EQUATIONS 

From  Equations  C.l  and  C.2,  the  informed  reader  may  see  that  these  ordinary 
differential  equations  can  be  solved  to  obtain  expressions  for  force  levels  as  a  function 
of  time.  Consider  the  following  differential  equations: 
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dx 

Tt   =  ~&y{t) 

i  -  -^ 

which,  when  differentiated  with  respect  to  the  independent  time  variable  t,  yield 


d2x  A  dy 

~dt2    ~  ~a~dt 

cPy  -dx 

~dP     =     ~  ~dt' 

dy  dx 

Substituting  the  values  for  —  and  — -  and  using  the  initial  conditions  asso- 

dt  dt 

ciated  with  Equations  C.l  and  C.2,  one  has  the  following  two  (second  order,  linear, 
constant  coefficient,  and  homogeneous)  ordinary  differential  equations  (ODE).  For  a 
more  detailed  explanation  of  ODEs  of  this  type,  consult  [Ref.  13]. 


d2x 

dt2 

= 

abx(t) 

x(0) 

= 

x0 

initial  condition  1 

dx(0) 
dt 

= 

-ay0 

initial  condition  2 

cPy 
dt2 

= 

bay(t) 

y(o) 

= 

Vo 

initial  condition  1 

4,(0) 

—  hx.n 

initial  condition  9. 

(C.3) 


(C.4) 


dt 

As  such,  each  second  order  ODE  can  be  solved  using  a  characteristic  equation 
after  assuming  general  exponential  solutions  of  the  form  x(t)  =  C\  expnt  +C2  expr2t 
and  y(t)  —  c3expr3< +c4expr4<.  Considering  the  second  order  ODE  Equation  C.3  and 
its  initial  conditions  in  isolation,  suppose  that  x(t)  =  exprt,  where  r  is  a  parameter 
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to  be  determined.  Then  it  follows  that  —  =  rexprt  and  -r-r  =  r2exprt.  Substituting 

at  at1 

these  values  into  Equation  C.3  yields 


2  rt  *  l  rt 

rexp       =    aoexp 
{r2-ab)exprt    =    0. 

And,  since  exprt  >  0  for  all  values  of  t,  then  (r2  —  ab)  must  equal  zero.  Therefore, 
r  =  ±yd&  are  the  roots  to  the  characteristic  equation.  Note  that  a  and  b  are  both 
positive  real  numbers  implying  that  ±y  ab  are  both  positive  real  numbers.  Knowing 
this,  the  complex  and  repeated  root  cases  can  be  eliminated  [Ref.  13].  General 
solutions  to  the  second  order  ODE  (Equation  C.3)  can  now  be  represented  by 

x{t)    =    c1exp(^f+c2exp-(V^)t. 
Examining  the  initial  condition  1, 

x(0)     =    x0 

<=>  C\  +  C2      =      Xo, 

and  initial  condition  2, 


dx(0) 

^r  =  ~ayo 

C\\J ab  —  C2\  ab    =     —ayo. 

The  reader  will  recognize  that  the  resulting  two  equations  representing  the 
initial  conditions  applied  to  our  general  solution  form  a  system  of  two  linear  equa- 
tions in  two  unknowns,  namely,  cx  and  c2.    Solving  this  system  of  equations  yields 

C\  —  —  I  Xq  —  \  -yo  I  and  c2  =  -  \  xo  -\-  \  -yo  ]  ■  Substituting  these  values  back  into 
2\  \  b  2\  V6 
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the  general  form  of  the  solutions,  Equations  C.5  and  C.6  result.  (Without  loss  of  gen- 
erality, Equation  C.6  follows  from  the  work  presented  for  Equation  C.5;  the  interested 
reader  is  encouraged  to  verify  the  result.) 


x(t)     =     ULo-JLAexp^+lxo  +  Jjyo 


V(t)    =    \  j  L  "  \j\xo\  exp^X  +  L  +  \j\x0  |  exp"^ 


(C.6) 


As  an  alternative  to  the  exponential  representation  of  the  force  size  equations 


(dependent  on  timet),  one  can  make  use  of  the  hyperbolic  functions  cosh  = 

exp   —  exp 
and  sinh  = to  transform  Equations  C.5  and  C.6  into 


exp  -(-  exp 


-6 


x(t)    =    x0  cosh 
y(t)    =    y0  cosh 


abjt 
ab]t 


—  vow  -  sinh 

lb 

ft  ■   u 

—  xq\\  -  sinn 

V  a 


abjt 
ab]t 


(C.7) 
(C.8) 


The  most  significant  advantage  to  this  form  is  the  insight  provided  for  the  parametric 
relationships  that  exist  within  the  force  level  equations  themselves.  Consider  the 
comparison  of  force  size  at  time  t  compared  to  the  initial  force  size  for  a  hypothetical 
force  X.  This  can  be  accomplished  by  dividing  Equation  C.7  by  xq,  resulting  in 


Xq 


=    cosh 


ab)  t 


yo    a 

—  \  —  sinh 
xo  V  b 


ab)t 


(C.9) 


From  Equation  C.9,  the  reader  can  see  that  the  level  of  force  X  appears  to  be  de- 
pendent upon  three  terms;  the  initial  force  ratio  (  —  ) ,  the  relative  fire  effectiveness 

\x0J 

measure  I  —  J ,  and  what  can  be  called  the  intensity  of  combat  indicator  (  y  ab  j .  In 
a  battle  to  annihilation,  these  indicator  values  can  be  used  to  generalize  simulation 
results  [Ref.  9]. 
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Hence,  one  can  now  deterministically  solve  for  force  sizes,  x(t)  and  y(t),  for 
the  desired  time  t  using  Equations  C.5  and  C.6  (or  Equations  C.7  and  C.8)  and  the 
estimated  attrition  rate  coefficients  a  and  b.  This  practical  extension  of  Lanchester's 
Square  Law  (F|F)  affords  the  analyst  a  straightforward  method  for  estimating  force 
size  throughout  a  battle.  The  force  size  results  can  also  be  used  to  verify  a  and  6 
themselves  by  comparing  the  numerical  values  of  x(t)  and  y(t)  to  the  analogous 
high  resolution  results  (force  size)  at  various  points  throughout  the  simulation.  This 
practical  exercise  is  what  remains  to  be  examined  in  this  Appendix. 

4.  COMPARISON  OF  JANUS  SIMULATION  RESULTS 
TO  THE  DETERMINISTIC  LANCHESTER  SQUARE 
LAW  RESULTS 

Using  the  values  for  a  and  b  derived  from  the  Janus  simulation  data  output, 
one  can  form  the  two  Lanchester  Square  Law  (F|F)  Equations  CIO  and  C.ll  below. 


dx 

—    =    -(0.00217256)y(*)       with  x0  =  113  (C.10) 

at 

J.    =    -(0.00072153)x(t)       withy0  =  55  (C.ll) 

at 

Or,  equivalently,  using  Equations  C.5  and  C.6  the  force  size  equations  can  be  written 
as 


x(t)    =    H(ll3-^55Jexp^)*+(ll3  +  jf55)exp-(^)*}    (C.12) 


where 


y(t)    =    -{\55-\l-m)exp^^t+\55  +  \l^m\exp-^^t\    (C.13) 


a    =    0.00217256 
6  =  0.00072153. 
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Calculating  the  force  size  values  of  x(t)  and  y(i)  for  values  of  t  =  (0, 15,30,  •  •  •  ,210) 
minutes  using  Equations  C.12  and  C.13  above,  the  reader  will  note  that  the  deter- 
ministic approximations  depicted  in  Table  XV  on  page  105  of  this  Appendix  are  quite 
representative  of  the  Janus  simulation  results.  (Note:  The  deterministic  results  for 
x(t)  and  y(t)  are  shown  as  integer  results  which  were  derived  by  rounding  the  contin- 
uous solutions  to  Equations  C.12  and  C.13).  Employing  the  deterministic  Lanchester 
equations  to  replicate  the  high  resolution  simulation  results  is  called  "replaying"  the 
simulation.  This  is  the  role  that  the  ATCAL  Phase  II  algorithm  plays  as  part  of  the 
aggregated  models  it  supports.  In  this  example,  the  deterministic  Lanchester  results 
estimate  the  casualties  achieved  during  the  Janus  simulation  with  fairly  high  accuracy 
as  depicted  by  Figures  11  and  12  on  page  106.  Note  that  the  maximum  absolute  error 
of  the  approximations  is  no  greater  than  two  weapon  systems  for  both  forces. 

The  small  discrepancies  can  be  attributed  to  several  facts.  First,  the  author  ar- 
bitrarily decided  to  single  out  two  opposing  weapon  systems  from  a  highly  integrated 
heterogeneous  combat  simulation  to  demonstrate  the  mathematics  associated  with  the 
homogeneous  form  of  Lanchester's  Square  Law  (F|F)  and  the  COMAN  MLE  attrition 
coefficient  estimation  process.  The  reader  may  have  already  gathered  that  the  high 
resolution  simulation  presented  in  Chapter  II  infers  the  need  for  a  heterogeneous-type 
attrition  analysis  construct;  however,  tractable  formulations  do  not  exist  for  het- 
erogeneous Lanchester-type  equations.  So,  in  presenting  the  simplistic  closed  form 
solutions  associated  with  the  homogeneous  problem,  the  reader  is  able  to  gain  an 
appreciation  for  the  mathematics  behind  attrition  computations.  Second,  examin- 
ing the  casualty  times  listed  in  Table  XIV,  one  can  see  that  the  exchanges  of  aimed 
fire  between  the  REDFOR  BMP-2  (with  AT-5)  and  the  BLUEFOR  M2  (with  TOW) 
were  focused  around  several  small  battles.  Each  force  shared  losses  during  several 
distinct  time  periods,  i.e.,  intervals  where  [tk  —  tk-i]  are  small  and  follow  each  other 
closely.  The  larger  values  of  [tk  —  tk-i]  depict  periods  of  inactivity  during  the  battle. 
Perhaps  distinct  COMAN  MLE  attrition  rate  estimates  a  and  b  could  be  derived  for 


102 


several  time  sub-intervals  representing  their  corresponding  sub-battles.  Lastly,  the 
reader  should  recall  that  this  research  only  utilizes  one  stochastic  sample  of  the  high 
resolution  simulation.  A  more  statistically  robust  analysis  could  have  been  made 
had  the  attrition  rate  parameters  been  estimated  after  considering  the  time  series 
Markov  (Poisson)  attrition  processes  for  several  independent  simulations. 

Even  with  the  perceived  shortcomings,  the  reader  will  agree  that  the  sampled 
deterministic  Lanchester  replay  of  the  high  resolution  simulation  is  quite  good.  Ad- 
ditionally, through  the  application  and  understanding  of  this  simplified  example,  the 
reader  should  have  a  thorough  appreciation  for  the  expanse  of  this  subject  matter. 
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Casualty 

Target       No.      Firer 

k 

Time  Interval 

No.  of  M2 

No.  of  BMP-2 

Time 

Lost 

tk  —  h-i 

(rrik-i) 

(ifc-i) 

16.41 

M2 

L        BMP-2 

1 

16:41 

113 

55 

16.56 

BMP-2 

L        M2 

2 

0:15 

112 

55 

17.42 

M2 

L        BMP-2 

3 

0:46 

112 

54 

18.27 

M2 

L        BMP-2 

4 

0:45 

111 

54 

18.28 

BMP-2 

L        M2 

5 

0:01 

110 

54 

37.26 

M2 

L        BMP-2 

6 

18:58 

110 

53 

37.49 

M2 

L        BMP-2 

7 

0:23 

109 

53 

38.10 

BMP-2 

L        M2 

8 

0:21 

108 

53 

38.11 

M2            : 

L        BMP-2 

9 

0:01 

108 

52 

46.32 

M2            : 

L       BMP-2 

10 

8:21 

107 

52 

49.24 

BMP-2 

L        M2 

11 

2:52 

106 

52 

59.29 

BMP-2       : 

L        M2 

12 

10:05 

106 

51 

64.14 

BMP-2 

[        M2 

13 

4:45 

106 

50 

83.12 

M2              1 

[        BMP-2 

14 

18:58 

106 

49 

84.20 

BMP-2       ] 

L        M2 

15 

1:08 

105 

49 

84.21 

M2              ] 

L        BMP-2 

16 

0:01 

105 

48 

110.58 

M2              ] 

L        BMP-2 

17 

2:37 

104 

48 

111.13 

M2              ] 

I       BMP-2 

18 

0:15 

103 

48 

131.51 

M2              ] 

[        BMP-2 

19 

20:38 

102 

48 

142.28 

BMP-2       ] 

[        M2 

20 

10:37 

101 

48 

142.30 

M2             ] 

L        BMP-2 

21 

0:02 

101 

47 

143.03 

BMP-2       ] 

[       M2 

22 

0:33 

100 

47 

143.10 

BMP-2       ] 

[       M2 

23 

0:07 

100 

46 

149.26 

M2              ] 

L       BMP-2 

24 

6:16 

100 

45 

149.30 

BMP-2       1 

L        M2 

25 

0:04 

99 

45 

149.48 

M2              ] 

L       BMP-2 

26 

0:18 

99 

44 

153.40 

M2              ] 

L        BMP-2 

27 

3:52 

98 

44 

173.36 

M2              ] 

L        BMP-2 

28 

19:56 

97 

44 

183.35 

BMP-2       1 

M2 

29 

9:59 

96 

44 

[Ref.  4] 


Table  XIV.  Janus  Simulation  Data:  Coroner's  Report 
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t 

Janus  x(t) 

Lanchester  x(t) 

Janus  y(t) 

Lanchester  y(t) 

0 

113 

113 

55 

55 

15 

111 

113 

54 

55 

30 

109 

110 

53 

53 

45 

108 

107 

51 

52 

60 

106 

106 

50 

50 

75 

105 

106 

49 

49 

90 

103 

104 

48 

48 

105 

101 

104 

47 

48 

120 

100 

102 

46 

48 

135 

98 

101 

45 

48 

150 

97 

98 

44 

44 

165 

96 

97 

43 

44 

180 

94 

96 

42 

44 

195 

93 

96 

41 

43 

210 

92 

96 

40 

43 

Table  XV.  Comparison  of  Janus  Attrition  Data  to  Lanchester- Type  Replay 
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Figure  11.  Graph  of  Force  Size  Over  Time  for  BLUEFOR  Simulation  vs.  Lanchester 
Replay  Results 
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Figure  12.  Graph  of  Force  Size  Over  Time  for  REDFOR  Simulation  vs.  Lanchester 
Replay  Results 
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